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I  .  INTRODUCTION 

Many  of  the  properties  of  the  ionospheric  polar  cap  and  auroral  zone  are  stnMiely 
affected  by  particle  and  energy  flow  between  the  ionosphere  and  the  magnetosphere,  in 
order  to  study  such  coupling  processes,  we  have  applied  both  kinetic  descriptions  of 
particle  transport  and  acceleration  and  fluid  descriptions  of  ambient  ionospheric  ion 
densities.  To  understand  and  model  the  physics  of  kinetic  phenomena,  we  have  used 
plasma  and  Monte  Carlo  simulations  as  well  as  numerical  transport  calculations.  Our 
approach  has  been  to  develop  these  models,  test  them  against  observations  and  then  refine 
them  accordingly.  But  often,  the  available  observations  are  not  of  the  quantities  calculated 
in  kinetic  models.  Rather,  it  is  the  effects  of  kinetic  phenomena  that  are  often  measured 
Thus,  in  order  to  test  and  refine  the  kinetic  models,  complementary'  models  are  needed  that 
will  calculate  how  the  kinetic  phenomena  affect  such  quantities  as  ion  densities,  ion  and 
electron  temperatures,  and  optical  emissions. 

In  sections  two  and  three,  we  describe  our  kinetic  modeling  of  the  flow  of  electrons 
between  the  ionosphere  and  magnetosphere.  The  next  two  sections  focus  on  coupling 
through  ion  transport  and  acceleration.  Our  fluid  modeling  of  the  E  and  F  regions  are 
treated  in  sections  six  through  eight.  A  final  section  summarizes  the  overall 
accomplishments  of  this  contract. 


2.  IONOSPHERIC  ELECTRON  FLUXES  IN  THE  DIFFUSE  AURORA 

The  nighttime  high  latitude  ionosphere  (auroral  ionosphere)  is  produced  by 
precipitating  particles  (electrons  and  ions)  with  kinetic  energies  ranging  from  .  1  to  more 
than  10  keV.  The  diffuse  aurora  is  a  rather  uniform,  relatively  low  intensity,  oval-shaped 
band  of  auroral  emissions  present  at  essentially  all  local  times,  displaced  equatorward  of, 
but  usually  overlapping,  the  oval  of  discrete  auroral  arcs.  The  diffuse  aurora  is  generated 
by  kilovolt  electrons  and  protons  precipitating  from  the  plasma  sheet  and  interacting  with 
the  neutral  atmosphere.  It  was  recognized  in  the  late  19^s  (Kennel,  1969)  that  the  spatial 
structure  of  the  diffuse  aurora  could  be  explained  as  a  consequence  of  the  coupling  of 
plasma  sheet  convection  towards  the  Earth  and  the  loss  of  trapped  particle  flux  by 
precipitation  into  the  ionosphere.  As  the  flow  nears  the  Earth,  convection  slows  while  the 
widening  loss  cone  and  falling  bounce  time  make  precipitation  more  effective.  To  analyze 
this  situation,  Kennel  introduced  a  powerful  approximation  to  describe  the  effect  of  pitch 
angle  scattering,  namely,  the  approximation  of  strong  diffusion,  that  states  that  an  isotropic 
pitch  angle  distribution  is  maintained  in  the  face  of  all  the  processes  that  tend  to  produce 
anisotropy.  This  approximation  has  wide,  but  not  complete,  observational  justification, 
and  is  of  such  power  in  simplifying  the  calculations  of  precipitation,  that  it  has  been  almost 
universally  adopted  in  calculations  of  this  kind  (Fontaine  and  Blanc,  1983;  Wolf,  1983). 
Examination  of  the  observations  in  more  detail,  however,  shows  that  electron  flux  in  the 
diffuse  aurora  is  not  always  isotropic  over  the  loss  cone  (Sharber,  1981)  and  that  the 
conditions  of  strong  pitch-angle  scattering  do  not  apply. 

We  have  developed  a  model  for  the  structure  of  the  diffuse  aurora  and  the 
precipitation  therein  using  a  full  kinetic  description  of  the  particle  fluxes,  so  that  the 
inaccurate  assumption  of  isotropy  can  be  abandoned.  By  including  the  effects  of 
convection,  precipitation,  and  backscatter  along  with  pitch  angle  scattering,  the  structure  of 
the  auroral  ionosphere  can  be  calculated  more  completely  and  with  a  greater  degree  of 
realism  than  before. 
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2.1  One  Dimensional  Kinetic  Model 


The  fundamental  equation  defining  the  problem  is  the  bounce- averaged  pitch-angie 
diffusion  equation  with  an  advective  term  describing  the  effects  of  plasma-sheet 
convection.  It  was  realized  that  eigenfunction  treatments  of  the  equation  would  not  be 
adequate,  and  that  an  initial-value  solution  was  necessary,  starling  with  a  plasma  velocity 
distribution  at  a  high  L-shell  and  letting  the  plasma  convect  Earthward.  We  have  developed 
the  algorithm  for  this  solution.  To  combine  the  effects  of  convection  and  diffusion  on  the 
velocity  distribution  function  over  a  time  step,  an  "operator-splitting"  technique  was  used  in 
which  separate,  stable  algorithms  for  advancing  the  velocity  distribution  function  by- 
convection  alone  and  by  diffusion  alone  are  alternately  applied  to  describe  their  combined 
effects.  The  convective  algorithm  is  an  integration-along-characteristic-trajectory  solution 
of  the  Liouville  equation,  assuming  that  an  electron's  magnetic  moment  and  longitudinal 
invariant  are  conserved  during  the  time  step.  The  guiding  center  drift  is  assumed  to  be 
E  X  B  drift.  Hence,  the  model  will  simulate  the  latitudinal  dependence  of  the  electron  flux. 
The  method  was  demonstrated  to  conserve  flux-tube  particle  content  over  finite  time  steps, 
so  long  as  the  energy  and  pitch-angle  grids  were  carefully  chosen;  in  particular,  the  energy 
grid  must  extend  to  high  enough  energies  to  allow  for  the  acceleration  due  to  "adiabatic” 
compression.  To  solve  the  diffusion  equation,  a  finite-difference  approximation  to  the 
pitch-angle  diffusion  term  was  worked  out,  with  special  care  taken  to  make  the  algorithm 
time-centered  to  help  insure  stability.  The  result  is  of  the  Crank- Nicholson  form.  The 
boundary  condition  on  the  pitch-angle  derivative  of  the  velocity  distribution  at  small  pitch 
angle  is  important  because  it  determines  the  flux  out  of  the  trapped  distribution  that  is  the 
net  precipitative  flux  into  the  ionosphere.  This  flux  is  determined  by  the  boundary-layer 
treatment  of  the  loss  cone  in  velocity  space,  coupled  with  a  transport  treatment  of 
atmospheric  backscatter  and  secondai7  production.  To  achieve  time-centering,  this  flux  is 
calculated  at  the  half  time  steps.  All  finite-difference  expressions  are  carried  out 
consistently  to  second  order  in  pitch  angle  spacing;  this  promotes  the  cancellation  of 
truncation  error  in  evaluating  integral  fluxes. 

When  pitch  angle  scattering  is  not  strong,  the  results  of  the  model  show  that 
precipitation  does  not  empty  a  convecting  flux  tube  before  it  has  drifted  past  the 
equatorward  edge  of  the  diffuse  aurora.  Thus,  some  other  mechanism  must  be  introduced 
to  explain  the  observed  drop-off  of  precipitation  at  the  edge.  Others  have  determined  that 
the  Earthward  boundary  of  the  convective  motion  of  particles  injected  from  the  tail  can 
explain  the  energy  dependence  of  the  boundary  in  a  number  of  events.  We  have 
incorporated  drift-boundary  limitation  of  flux  into  the  kinetic  model.  One  benefit  of  having 
the  equatorward  edge  of  the  diffuse  aurora  limited  only  by  drift  boundaries  is  that  in 
individual  events,  observation  of  the  boundary  then  gives  us  a  direct  handle  on  the 
magnitude  of  the  convection  electric  field;  when  the  boundary-  is  detemiined  by  exhaustion 
of  the  flux  tube,  the  position  of  the  boundary  depends  on  the  combination  of  the  effects  of 
the  strength  of  pitch  angle  scattering  and  the  drift  velocity,  and  no  direct  means  of 
determining  the  convection  field  through  the  precipitation  pattern  exists. 


2.2  Two  Dimensional  Kinetic  Model 

The  two  dimensional  version  of  our  model  requires  that  the  guiding  center  drift  be  a 
combination  of  E  X  B  and  gradient/curvature  drifts.  This  allows  the  mcxlel  to  map  out  the 
entire  two  dimensional  structure  of  the  band  of  diffuse  aurora  precipitation.  At  present,  the 
model  is  too  computationally  expensive  to  be  applied  to  the  large  number  of  trajectories 
needed  to  resolve  the  details  of  both  the  invariant-latitude  and  local-time  variations  of  the 
flux,  and  only  the  coarsest  local-time  variations  have  been  admitted  into  the  model.  This 
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has  been  justified  by  the  observed  smooth  local-time  structure  of  the  diffuse-auroral 
precipitation. 

Nevertheless,  statistical  studies  of  the  precipitation  do  find  significant  local-time 
structure  which  it  would  be  interesting  and  important  to  model  in  full  detail.  As  a  first  step 
toward  such  a  model,  we  have  developed  a  test-bed  two  dimensional  plasma-sheet 
convection  model  with  which  to  explore  the  consequences  of  the  two  dimensional  nature  of 
plasma  sheet  convection.  This  model  uses  the  popular  assumption  of  isotropic  pitch -angle 
distributions,  but  is  structured  in  a  way  similar  to  our  full  kinetic  model  so  that  as  more 
efficient  algorithms  for  the  pitch-angle  scattering  and  atmospheric  tran.sp>on  are  developed, 
they  can  be  incorporated  into  the  two-d  model  for  realism.  As  a  first  application  of  the  tv-o- 
d  model,  a  loss-less  flow  was  studied,  demonstrating  how  the  limits  on  inward  convection 
on  the  nightside,  caused  by  the  energy-dependent  gradient/curvature  drifts  of  the  electrons, 
can  produce  a  realistic  looking  equatorward  edge  for  the  diffuse  aurora. 


2.3  Future  Directions 

With  a  kinetic  model  for  the  electron  flux  in  the  diffuse  aurora  developed,  several 
applications  of  the  model  can  be  considered.  Possible  future  investigations  include  the 
following;  considerable  insight  into  the  form  of  the  electron  flux  can  be  obtained  by 
studying  the  ways  in  which  it  varies  when  the  functional  inputs  to  the  model,  such  as  the 
initial  flux,  the  diffusion  coefficient,  and  the  electric  potential,  are  varied.  Such  a 
systematic  study  of  these  variations  is  obviously  required.  Despite  the  complexity  of  the 
phenomena  that  control  the  electron  flux,  the  resulting  form  is  of  remarkable  simplicity. 
From  the  complicated  output  of  the  model,  we  would  like  to  be  able  to  extract  a  simple 
analytical  explanation  of  the  electron-flux  form.  Such  results  could  be  incorporated  into 
global  models  of  the  inner  magnetosphere.  Finally,  because  one  portion  of  the  model 
de.scribes  the  atmospheric  transport  of  the  precipitating  electrons,  it  is  to  be  noted  that 
among  the  results  of  the  model  is  the  energy  deposited  in  the  ionosphere  by  the 
precipitating  electrons.  This  quantity  can  be  used  to  estimate  ionospheric  conductivities, 
optical  emissions,  and  the  electron  density  profile  in  the  E  region  of  the  ionosphere. 


3.  ENERGETIC  PHOTOELECTRONS  AND  POLAR  CAP 

PRECIPITATION 

The  Low  Altitude  Plasma  Instrument  (LAPI)  onboard  the  Dynamics  Explorer  2 
satellite  has  observed  an  energetic  electron  population  (with  energies  between  60  and  850 
eV)  in  the  daytime  topside  ionosphere.  This  energetic  population  exists,  together  with  the 
"classical"  photoelectrons  (at  energies  less  than  ^  eV)  observed  by  earlier  satellites  and 
rocket  experiments.  Preliminary  analysis  of  this  data  suggests  that  this  energetic  population 
is  produced  by  soft  solar  X-ray  fluxes  in  the  same  way  that  EU  V  (extreme  ultraviolet)  solar 
fluxes  produce  the  classical  photoelectrons  (Winningham  et  al.,  1989).  Further  work 
indicates  that  both  the  classical  and  the  energetic  photoelectrons  are  present  in  the  sunlit 
polar  cap.  The  focus  of  our  work  is  to  model  the  electron  flux  in  the  polar  cap  and 
discriminate  between  this  photoelectron  population  and  the  well-known  polar-rain 
population  of  magnetospheric  origin.  However,  before  modeling  polar  cap  data  in  detail, 
we  must  return  to  the  midlatitude  data  in  order  to  settle  some  unresolved  issues  from  the 
preliminary  analysis. 
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з. i  Midlatitude  Calculation 

In  modeling  midlatitude  electron  flux  measurements  by  photoelectron  transport 
calculations,  the  major  issue  is  the  uncertainty  in  the  spectrum  of  the  solar  EU  V  and  soft  X- 
ray  flux.  The  solar  flux  is  one  of  the  primary  input  parameters  of  such  models  because  its 
magnitude  and  spectrum  determines  the  magnitude  and  variation  with  energy  of  the 
photoeleclTon  source  function.  Because  of  the  temporal  variability  of  the  solar  EUV  and 
soft  X-ray  spectrum,  considerable  uncertainty  might  exist  in  an  ab-initio  calculation  of  the 
electron  flux  in  any  individual  event.  We  resolved  this  issue  by  turning  the  question 
around:  by  requiring  that  modeled  electron  fluxes  agree  with  the  satellite  observations,  we 
could  deduce  a  possible  form  for  the  solar  spectrum.  This  synthetic  spectrum  could  then  be 
used  to  model  the  photoelectrons  in  electron  data  taken  in  regions  where  photoelectrons 
were  not  the  only  constituents,  such  as  the  sunlit  polar  cap.  In  general,  this  was  a 
successful  approach,  in  the  sense  that  over  most  of  the  wavelength  range,  the  synthetic- 
spectrum  did  not  depan  too  far  from  the  statistical  expectations  for  the  spectrum  under  the 
lesel  of  solar  activity  found  during  the  event  modeled.  In  order  to  obtain  detailed 
agreement  with  elecnon  flux  data  between  300  and  600  eV,  however,  we  initially  had  to 

и. se  fairly  large  solar  soft  X-ray  fluxes,  especially  between  20  and  40  Angstroms  (A). 
Improvements  and  refinements  made  to  the  neutral  atmosphere  model,  the  photoionization 
cross  sections,  the  photoabsorption  cross  sections,  and  the  model  parameters,  only  made 
the  problem  worse.  Although  the  solar  X-ray  spectrum  is  highly  variable,  it  was  not  at  all 
c'ear  that  our  inferred  soft  X-ray  spectrum  between  20  and  40  A  was  physically  rettsonable. 

This  difficulty  prompted  a  search  for  processes  that  might  produce  "more"  high 
energy  electrons.  One  process  we  considered  was  the  production  of  Auger  electrons.  We 
found  that  for  O,  N2,  and  O2,  a  K-shell  ionization  is  followed  over  99%  of  the  time  by  an 
.Auger  transition  that  produces  another  free  electron.  Further,  the  energy  of  that  electron  is 
dependent  on  the  neutral  species  involved  but  independent  of  the  photon  energy.  Thus, 
any  photon  with  enough  energy  to  cause  a  K-sheli  ionization  will  lead  to  the  production  of 
a  second  electron,  the  Auger  electron,  within  a  narrow  energy  band  at  high  energy.  When 
we  included  this  process  in  our  model,  we  found  two  "bumps"  appearing  in  the 
photoelectron  flux  exactly  where  similar  features  are  seen  in  the  electron  flux  data.  The  soft 
solar  X-ray  flux  needed  to  fit  the  data  was  significantly  smaller  than  that  needed  when 
Auger  electrons  were  not  included  in  the  calculation. 


3.2  Polar  Cap  Calculation 

Having  greater  confidence  in  our  midlatitude  results  and  the  deduced  synthetic  solar 
spectrum,  we  used  that  solar  spectrum  in  a  series  of  calculations  designed  to  evaluate  the 
role  of  the  photoelectrons  in  the  sunlit  polar  cap.  In  these  calculations,  we  considered  two 
different  levels  of  flux  fur  the  precipitating  polar  rain  electrons:  a  typical  or  moderate  level, 
as  characterized  in  the  statistics  of  the  magnitudes  of  flux  observed  by  DMSP,  and  a  weak 
level  representing  the  minimal  polar  rain.  For  each  case,  we  performed  three  calculations: 
a  full  calculation  with  both  photoelectron  and  polar  rain  sources  of  electrons,  a  calculation 
with  only  photoelectron  production  fno  polar  rain),  and  a  calculation  with  Just  polar  rain 
(no  photoproduction).  For  both  minimum  and  moderate  polar  rain,  we  found  that  there  can 
easily  be  conditions  where  upflowing  electrons  in  the  polar  cap  are  dominated  by 
photoeiecirons  at  all  altitudes.  For  downflowing  electrons,  on  the  other  hand,  the 
importance  of  the  photoelectrons  is  a  very  .sensitive  function  of  altitude.  At  8(X)  km.  it 
appears  that  the  downflowing  electrons  can  be  treated  purely  as  a  polar  rain  population. 
When  the  polar  rain  is  minimal,  the  photoeleclron  population  at  a  representative  energy  of 
100  eV  can  become  important  at  600  km  while  for  moderate  polar  rain,  this  transition  can 
occur  around  300  km.  These  transition  altitudes  also  depend  on  the  energy  being 
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considered.  What  does  this  mean  for  polar  rain  studies  in  the  sunlit  polar  cap?  if  the  data 
come  from  a  zenith  looking  insn-ument  onboard  a  satellite  at  HIM)  km,  such  as  DMSF.  theti 
the  "contamination"  by  ionospheric  photoelectrons  will  be  minimal  and  the  data  can  be 
considered  "pure”  polar  rain.  On  the  other  hand,  data  from  a  satellite  such  as  DE-2,  which 
went  as  low  in  altitude  as  300  km,  may  very  well  include  a  significant  photoeleciron 
population,  and  the  relative  contributions  of  polar  rain  precipitation  and  photoelectron 
production  would  have  to  be  considered  carefully  on  a  case-by-case  basis. 

To  this  point,  our  modeling  of  ionospheric  photoelectron  flux  data  had  been  done 
on  a  case-by-case  basis,  that  is  to  say,  data  from  different  locations  had  been  mcxlclcd 
separately. 

We  generalized  our  modeling  of  low-altitude  electron  fluxes  at  moderate  energies 
(,20  eV  to  1  keV)  to  be  able  to  calculate  the  fluxes  along  a  satellite  track  from  the  daytime 
midlatitudes,  through  the  cusp,  over  the  polar  cap,  and  across  the  nighttime  auroral  region. 
In  this  way,  we  were  able  to  simulate  satellite  observations  of  the  electron  flux  through  an 
entire  segment  of  an  orbital  pass  and  compare  the  theoretical  results  to  the  satellite 
observations.  Using  this  global  model,  we  simulated  electron  flux  data  taken  on  a  DE-2 
overflight  of  the  northern  polar  cap.  The  electron  data  came  from  the  Low  Altitude  Plasma 
Instrument  (LAPl)  onboard  DE-2,  from  the  same  pass  that  was  used  in  our  earlier  work  on 
energetic  photoelectrons  in  the  midlatitudes.  For  this  simulation,  the  solar  flux  inferred  by 
modeling  electron  data  at  one  midlatitude  location  was  used  in  calculating  the  photoelectron 
source  function  throughout  the  pass,  while  the  satellite  observations  of  downfiowing 
electron  fluxes  were  used  as  a  boundary  condition  on  the  region  in  which  we  were  solving 
the  electron  transport  equation.  The  accuracy  of  the  modeling  can  then  be  evaluated  by 
comparing  the  calculated  upflowing  electron  flux  as  a  function  of  pitch  angle  and  energy  to 
the  flux  observed  by  satellite.  At  high  energies  (85-640  eV),  we  obtained  excellent 
agreement  between  data  and  theory  at  all  locations  along  the  satellite  track,  including  the 
terminator  region.  At  low  energies  (20-85  eV),  we  obtained  qualitative  agreement  with  the 
shape  of  the  latitudinal  electron  flux  profile,  but  we  underestimated  the  magnitude  of  the 
electron  flux.  Because  the  agreement  worsens  as  photoelectrons  make  up  a  smaller  fraction 
of  the  upflowing  flux,  we  suspected  that  our  treatment  of  low  energy  electron  backscatter 
needed  improvement.  In  particular,  uncertainties  in  backscatter  cross  sections  were  thought 
to  be  a  likely  source  of  these  errors. 

In  the  electron  transport  model,  the  essential  mechanism  for  angular  scattenng  of 
electrons  is  the  process  of  electron- neutral  elastic  collisions.  In  the  model's  formulation, 
the  differential  cross  section  for  electron-neutral  elastic  scattering  (as  a  function  of  energy 
and  pitch  angle)  is  taken  to  be  the  total  cross  section  times  a  normalized  phase  function, 
which  gives  the  probability  of  scattering  from  one  pitch  angle  to  another.  The  form  of  the 
phase  function  was  originally  derived  from  the  screened  Rutherford  formula  for  the 
differential  cross  section  (as  a  function  of  energy  and  scattering  angle).  Using  this  formula 
leads  to  an  angular  cross  .section  which  is  generally  forward  peaked.  An  examination  of 
experimental  cross  section  data  shows  that  this  is  a  reasonable  formulation  above  1(X)  eV 
but  the  cross  sections  below  100  eV  show  a  sizeable  additional  peak  in  the  backward 
direction.  For  example,  for  scattering  with  N2  at  energies  below  50  eV,  the  Rutherford 
formula  originally  u.sed  in  the  model  leads  to  a  probability  of  backward  scattering  that  is 
over  an  order  of  magnitude  smaller  than  that  seen  in  the  laboratory. 

To  examine  the  sensitivity  of  the  electron  transport  calculation  to  the  form  of  the 
phase  function,  we  made  a  test  calculation  using  an  isotropic  scattering  approximation 
below  100  eV,  This  is  still  a  poor  approximation  to  the  measured  cross  sections  but  we 
expected  it  to  give  a  significantly  larger  backscatter  compared  to  our  original  phase 
function.  In  this  way,  we  could  get  some  idea  whether  improving  our  phase  function 
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would  be  a  useful  modification  to  our  model.  In  the  sunlit  polar  cap,  the  test  calculation 
showed  only  small  changes  from  the  previous  results  (as  expected,  because  the  low-encrg , 
ilux  here  w  due  to  phoiopnxiuction  and  inelastic  scattering).  Through  the  icrminatuf, 
agreement  between  model  and  observations  improved  markedly.  In  darkness  some 
improvement  was  seen,  but  the  model  electron  fia.t  was  still  significantly  below  the 
observation:'.  Our  conclusion  was  that  the  backscaticr  is  indeed  quite  sensitive  to  the  fonn 
of  the  phase  function  and  that  modifying  our  fonnulation  of  it  would  be  worth  the  effort. 

Accordingly,  we  modified  the  electron  transport  model  to  use  a  more  precise 
representation  of  the  angular  dependence  of  the  clastic  enrss  section.  We  chose  to  follow 
the  approach  of  Prasad  et  al.  ( 1985)  who  used  the  parameterired  elastic  cross  sections  ot 
Poner  and  Jump  (1978).  'fhe  analytic  cross  sections  consist  of  two  screened  Rutherford 
type  terms,  and  depend  on  four  energy-dependent  parameters  deiemiined  by  nonlinear  least 
squares  fitting  of  measured  differential  cross  sections.  The  form  of  the  cross  sections 
allows  an  integration  over  pitch  angle  to  be  done  analytically  and  thus.  dtK*s  not  increase  the 
computer  requirements  of  the  model. 

The  electron  transport  mcxlel  is  actually  divided  into  two  codes.  The  first  is  called 
MX  and  it  prepares  the  cross  sections  into  a  matrix  form  for  use  in  solving  the  transport 
equatio'  .  The  second  s  named  B3C  and  it  actually  solves  the  transport  equation.  The 
changes  in  our  treatment  of  the  elastic  cross  sections  required  modifications  to  MX  hut 
required  no  changes  to  IV^C. 

When  comparison  with  data  was  made  using  the  new  version  of  the  model,  we 
found  that  at  low  energies  ('x:lcw  85  eV),  we  still  underestimated  the  upward  electron  fiux. 
f-'or  example,  at  iO  eV,  the  measurements  showed  an  electron  flux  of  6.0  x  lO-^  (el  enr^  s  ' 
sr’’  eV'l).  The  original  version  of  the  model  gave  a  flux  of  2.0  x  lO**  w'hile  the  tncxiificd 
model  gave  2.2  x  10-“’.  Clearly,  we  still  had  a  low  energy  problem.  We  then  made  a  series 
of  calculations  to  inve.stigate  the  role  of  several  "parameters"  in  the  mixlel,  specifically  the 
assumed  electron  density  and  the  assumed  pitan  angle  distribution  of  the  secondary 
electrons.  We  found  that  if  we  turned  off  the  energy  loss  due  to  Coulomb  scattering,  the 
electron  flux  at  10  eV  went  to  3.4  x  lO^  but  at  30  eV,  there  was  no  effect.  Thus,  a  different 
assumed  electron  density  would  do  little  to  solve  our  problem.  In  another  calculation,  we 
assumed  that  the  secondary  electrons  above  10  eV  were  all  forward  scaitered.  Tins  was  in 
contrast  to  our  normal  as.sumption  that  below  1(X)  eV,  the  secondaries  are  produced 
isotropically.  The  results  were  unexpected:  the  upgoing  flux  at  10  eV  increased  to  3.5  x 
lO^,  a  59%  increa.se  over  our  standard  calculation.  At  lower  Jtiiudes,  we  saw  an  increase 
of  over  a  factor  of  four  compared  to  our  standard  calculation.  The.se  results  were  so 
surprising,  we  decided  to  comptu'e  them  with  re.sults  from  other  versions  of  the  MX  and 
B3C  programs.  From  those  comparisons,  we  discovered  that  our  electron  transport  mtxiel 
would  give  dramatically  different  results  below  100  eV  when  different  pitch  angle  gnds 
were  used.  TTe  critical  difference  appeared  to  be  wheihi  r  the  grid  was  evenly  spaced  in  the 
cosine  of  the  pitch  angle  or  was  non-uniformly  spaced.  This  i,ssue  remains  unresolved,  but 
it  is  very  possible  that  when  I'csolved  the  problem  of  underestimating  the  low  energy  data 
w  ill  be  solved.  This  based  on  the  fact  that  our  previous  compariso  .s  to  data  involved 
calculations  on  a  non-unifonn  pitch  angle  grid,  whereas  when  we  use  the  evenly  spaced 
pitch  angle  grid,  we  calculate  a  flux  at  10  eV  of  5.0  x  10-.  This  is  well  within  the  error 
bars  of  the  data. 
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3.3  Future  Directions 


Our  work  on  electron  transport  could  lead  us  jn  a  variety  of  directions,  l-irst.  ihc 
type  of  analysis  that  we  performed  on  a  segment  of  one  DE-2  orbit  could  be  dune 
systematically  on  many  DE-2  orbits.  The  purpose  would  be  to  'est  the  aimo.sphera  - 
transport  methods  of  Strickland  et  al.  (1976)  under  a  variety  of  geophysical  condition'.. 
Another  topic  would  involve  coupling  our  iranspon  model  with  a  Monte  Carlo  simulation 
of  electron  scattenng  in  the  plasmasphere.  The  idea  would  be  to  develop  an  electron  kinetic 
model  for  an  entire  midlatitude  fluxtube.  The  role  of  the  Monte  Carlo  simulation  would  be 
to  simulate  adiabatic  motion,  Coulomb  collisions,  and  wave-panicle  interactions. 


4.  PROTON-HYDROGEN  ATOM  PRECIPITATION 

In  recent  years,  the  first  self-consistent  theory  for  the  combined  electron-proton - 
hydrogen  atom  aurora  was  developed  by  D.J.  Strickland,  B.  Basu,  J.R.  Jasperse.  and 
R.E.  Daniell.  The  theory  involves  the  use  of  transpon  theoretic  methods  to  solve  a  set  of 
three  Boltzmann  equations  for  the  elecucn,  proton,  and  hydrogen  atom  differential  fluxes. 
The  protons  and  hydrogen  atoms  are  coupled  to  each  other  by  charge-changing  collisions 
(charge  exchange  and  stripping),  and  the  electrons  are  coupled  to  the  protons  and  hydrogen 
atoms  because  of  the  secondary  electrons  that  are  produced.  The  new  feature  of  this  model 
over  previous  work  is  that  the  secondary  electrons  produced  by  the  primary  protons  and 
hydrogen  atoms  are  included  in  the  solution  of  the  electron  transport  equation.  The  major 
new  and  unexpected  result  from  doing  this  is  that  the  secondary  electron  flux  produced  by 
the  proton-hydrogen  atom  aurora  is  much  softer  than  that  produced  by  the  electron  aurora. 
This  is  due  to  the  fact  that  the  cross  sections  for  proton  and  hydrogen  impact  ionization 
decrease  exponentially  with  increasing  secondary  electron  energy,  whereas  the  cross 
sections  for  electron  impact  ionization  decrease  as  an  inverse  power  law  with  increasing 
secondary  energy. 


4.1  Comparisons  with  Observations 

Previously,  the  portions  of  the  model  involving  just  electron  transport  have  been 
tested  against  a  variety  of  data.  The  coupled  proton-hydrogen  atom  segment  of  the  model 
has  had  only  one  tvpe  of  comparison  with  data.  That  is,  a  calculation  of  the  ionization 
produced  by  a  specific  proton  aurora  agreed  with  the  ionization  that  was  observed  by 
incoherent  scatter  radar.  Until  the  work  reported  here,  the  coupling  between  the  electrons 
and  the  protons  and  hydrogen  atoms  was  untested. 

The  greate.st  difficulty  in  making  such  a  test  was  finding  a  .suitable  data  set.  Ideally, 
you  would  like  a  measurement  of  all  three  species  as  a  function  of  altitude,  energy,  and 
pitch  angle,  but  such  a  case  was  not  found.  Rather,  we  took  advantage  of  our  access  to, 
and  experience  with,  the  measurements  made  by  the  Low  Altitude  Plasma  Instrument 
(LAPI)  that  was  onboard  the  Dynamics  Explorer  2  satellite.  This  instrument  measured 
electron  and  ion  fluxes  as  a  function  of  energy  and  pitch  angle.  Such  data  can  serve  as  a 
limited  but  useful  test  of  the  kinetic  model.  In  particular,  it  allowed  us  to  test  how  well  the 
model  predicts  the  upgoing  electron  flux  that  results  from  incident  electrons  and  protons. 

As  stated  above,  our  primary  interest  was  in  the  e’ectrons  that  are  produced  by 
incident  protons.  Our  original  thought  was  to  look  at  three  cases;  the  pure  proton  aurora, 
the  pure  electron  aurora,  and  the  mixed  aurora.  But  after  studying  the  results  of  some 
preliminary  calculations,  we  came  to  realize  that  the  interesting  case  would  be  the  pure 
proton  aurora  as  compared  to  the  pure  electron  aurora.  The  mixed  aurora  is  not  very 
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uuercsiing  because  the  soft  electron  spectrum  produced  by  precipitating  protons  and 
hydrogen  atoms  can  easily  be  masked  by  the  electrons  produced  by  the  precipitating 
electrons.  This  also  meant  that  in  order  to  study  the  electrons  produced  by  precipitating 
protons  and  hydrogen  atoms,  we  would  need  a  very  "pure"  proton  aurora,  that  is, 
precipitating  protons  and  hydrogen  atoms  with  a  minimum  of  coincident  precipitating 
electrons.  This  need  for  a  very  pure  proton  aurora  made  finding  a  suitable  case  quite 
difficult.  For  example,  a  case  used  in  an  E  region  study  that  had  been  identified  previously 
as  a  pure  proton  aurora  did  not  come  close  to  being  pure  enough  for  studying  the  electroti 
backscatter  in  the  F  region.  We  realized  that  we  had  to  be  very  particular  where  we  kxiked 
for  proton  auroras.  Besides  looking  in  the  dusk  sector  equatorward  of  electron  auroras,  we 
also  focused  on  times  of  the  year  near  equinox.  This  was  done  in  order  to  find  a  ca.se 
w  here  the  ionospheric  portions  of  the  fieldline  would  be  in  darkness  so  that  there  would  be 
no  contamination  due  to  photoelectrons  produced  in  the  local  or  conjugate  ionosphere. 

We  found  a  suitable  case  from  day  289  in  1981  when  DE-2  passed  through  an 
electron  aurora  between  19:4.^  and  19:44  UT  and  then  passed  through  a  fairly  pure  proton 
aurora  between  19:44  and  19:45  UT.  The  solar  kx;al  time  was  around  21.3  hours  and  the 
satellite  altitude  was  around  950  km.  At  19:43,  the  satellite  was  at  an  invariant  latitude  of 
72  5  degrees  and  at  19:45,  it  was  at  66.5  degrees. 

This  case  on  day  289,  while  the  "best"  we  found,  was  far  from  ideal.  The  proton 
flux  was  quite  modest  requiring  that  we  sum  over  27  sweeps  (27  seconds)  to  obtain 
reasonable  counting  statistics.  Funher,  it  was  not  a  pure  proton  aurora  ca.se  but  rather  there 
was  a  measurable  flux  of  downgoing  electrons.  A  flux  of  upgoing  electrons  that  we 
wished  to  compare  to  our  model  was  ^so  measured  but  this  flux  was  also  fairly  modest  in 
intensity.  So  in  both  the  electron  and  ion  data,  we  were  forced  to  live  with  the  limitations 
of  low  counting  rates. 

Once  our  data  processing  was  in  place,  we  turned  to  modeling  our  case  from  day 
289.  As  noted  above,  the  measurements  were  made  around  950  km.  While  we  expected 
that  the  effects  of  charged  particle-neutral  collisions  would  be  small  at  this  high  altitude,  we 
wanted,  if  possible,  to  calculate  them  explicitly.  But  ..ttempts  to  execute  our  proton- 
hydrogen  atom  code  above  600  km  ran  into  numerical  problems.  After  some  work,  we 
came  to  believe  that  the  problem  was  not  an  error  as  such  but  rather  a  shoncoming  in  the 
design  of  the  algorithm  used  to  solve  the  coupled  proton  hydrogen  atom  transport 
equations.  Fortunately,  we  had  available  closed  form  analytic  solutions  for  this  problem 
that  are  appropriate  at  high  altitudes.  We  thus  used  the  analytic  solutions  from  950  km  to 
600  km  and  at  600  km,  these  solutions  were  used  as  the  boundary  conditions  for  a  full 
numerical  solution  at  altitudes  below  600  km. 

For  the  electron  transport  part  of  the  model,  we  decided  to  limit  ourselves  to 
calculations  from  600  km  down.  This  was  again  due  to  numerical  problems.  From  earlier 
work  we  knew  that  the  electron  transport  code  can  begin  to  lose  convergence  at  lower 
energies  as  you  attempt  to  solve  over  larger  ranges  of  altitudes.  But  based  on  our  previous 
experiences,  we  felt  it  was  useful  to  attempt  modeling  the  electron  data  as  if  it  were  taken  at 
6()()  km. 

From  our  first  comparisons,  we  discovered  that  the  model  predicts  upflowing 
electrons  that  are  certainly  the  right  order  of  magnitude.  Second,  we  found  that  the  model 
results  are  sensitive  to  the  proton  flux  above  27  keV.  But  LAPI  only  measured  ion  fluxes 
up  to  27  keV  and  thus,  our  results  are  sensitive  to  how  we  extrapolate  the  incident  proton 
flux  to  higher  energies.  Fitting  a  kappa  distribution  to  the  data  versus  a  Maxwellian  will 
lead  to  almost  a  factor  of  2  increase  in  the  inferred  incident  proton  fiiix.  There  is  evidence 
from  other  satellite  measurements  that  in  the  plasma  sheet,  a  kappa  distribution  provides  a 
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much  better  fit  to  data  than  does  a  Maxwellian  (Christen  et  al.,  19X8;  Christon  et  al.,  !98s): 
Christen  et  al.,  1991).  For  our  case,  using  a  kappa  extrapolation  gives  better  agreement 
with  data.  A  calculation  using  only  the  observed  incident  electron  flux  clearly 
underestimates  the  low  energy  data  (about  a  factor  of  three  at  50  eV),  Above  2(X)  eV.  the 
protons  have  little  effect  and  the  upgoing  electron  flux  results  from  the  incident  electrons. 
If  we  attempt  to  match  the  data  by  increasing  the  incident  electron  flux  above  50  eV,  wc 
have  to  use  a  flux  which  is  well  above  the  observed  incident  electron  flux.  I'hus.  the 
effects  of  the  protons  are  clearly  important  in  this  case.  Given  the  uncertainties  in  the  data 
and  the  need  to  extrapolate  the  proton  flux,  we  believe  we  can  claim  the  model  is  prixlucing 
results  good  to  a  factor  of  2. 

To  clearly  illustrate  the  contrast  between  upflowing  electrons  produced  by 
precipitating  protons  and  hydrogen  atoms  and  those  produced  by  precipitating  electrons, 
we  made  two  calculations.  In  one  case,  we  used  an  incident  flux  of  electrons  characterized 
by  a  1  keV  Maxwellian  with  a  total  energy  flux  of  1  erg  cm'^  sec  '.  For  the  other  ca.se,  we 
considered  an  incident  flux  of  protons  characterized  by  an  8  keV  Maxwellian  w  iih  a  total 
energy  flux  of  1  erg  cm‘2  secF  In  the  proton  case,  a  characteristic  energy  of  8  keV  was 
chosen  because  it  had  been  found  from  earlier  calculations  that  an  8  keV  proton  Maxwellian 
produces  a  volume  ionization  rate  which  peaks  at  about  the  same  altitude  as  that  produced 
by  a  1  keV  electron  Maxwellian.  When  we  compared  the  two  cases,  the  resultant  energy 
spectra  of  the  upflowing  electrons  were  dramatically  different.  Compared  to  the  electron 
case,  the  proton  case  produces  a  much  softer  spectrum,  a  spectrum  that  for  energies  above 
200  eV  is  over  4  orders  of  magnitude  smaller  than  in  the  electron  case.  Interestingly,  if  you 
examine  the  total  ionization  rates  for  the  two  cases,  you  see  that  they  differ  only  by  around 
25%. 


Using  data  from  the  Low  Altitude  Plasma  Instrument  (LAPI)  that  was  onboard  the 
DE-2  satellite,  we  compared  the  "purest"  proton  aurora  that  we  found  to  a  pure  electron 
aurora  that  had  about  the  same  incident  energy  flux.  We  found  clear  qualitative  agreement 
with  the  theoretical  prediction.  In  the  proton  case,  the  upflowing  electrons  above  KX)  eV 
were  clearly  a  couple  of  orders  of  magnitude  smaller  than  those  seen  in  the  electron  aurora 
case. 


One  step  in  improving  our  detailed  quantitative  modeling  of  data  from  day  289  in 
1981  was  to  extend  our  electron  transport  upper  boundary  from  600  km  to  900  km.  As 
expected,  this  led  to  convergence  problems  at  the  lowest  energies  considered,  but  we  were 
able  to  verify  that  above  30  eV,  collisions  have  very  little  effect  on  the  upgoing  electron 
flux.  Below  30  eV,  there  was  the  expected  attenuation  due  to  Coulomb  collisions  with  the 
thermal  electrons.  Thus,  detailed  comparisons  with  data  in  this  lower  energy  range  require 
both  an  accurate  electron  density  profile  as  well  as  an  accurate  treatment  of  the  Coulomb 
collisions. 

When  we  model  the  low  energy  upflowing  electrons  at  900  km,  where  do  they 
come  from?  Does  it  tell  us  anything  about  hovz  the  model  is  doing  in  the  E  region  where 
most  of  the  energy  deposition  takes  place?  To  address  these  questions,  a  set  of  calculations 
was  made  where  we  tried  different  lower  boundaries  from  98  to  2()0  km.  From  these 
calcjilations,  we  found  that  over  90%  of  the  upflowing  electrons  below  100  eV  come  from 
above  200  km.  Thus,  modeling  the  upflowing  electrons  seen  in  this  particular  event  on  day 
289  is  a  test  of  how  the  model  does  in  the  F  region  but  says  little  about  its  E  region 
performance. 
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4.2  Future  Direction.s 


Our  present  proton-hydrogen  atom  mtKicI  contains  no  processes  that  change  the 
pitch  angle  of  the  panicles.  Thus,  given  an  incident  downward  tlux.  of  protons,  the  mcxiel 
produces  no  upward  flux  of  protons  and  hydrogen  atoms.  On  the  other  hand,  an 
examination  of  ion  measurements  made  by  LAPl  during  a  proton  aurora  show  s  an  upward 
flux  of  ions  at  pitch  angles  near  90  degrees.  This  suggests  that  one  useful  improvement  to 
the  model  would  be  to  include  magnetic  mirroring  of  the  protons. 

Another  obvious  need  is  for  more  comparisons  with  observations.  In  this  work, 
we  have  examined  only  a  single  case  of  electrons  produced  by  precipitating  protons. 
Further,  there  have  been  no  comparisons  between  observations  and  predicted  proton  and 
hydrogen  atom  fluxes.  For  any  such  studies,  we  suspect  that  the  greatest  problem  will  ix- 
finding  suitable  data.  Finally,  the  optical  emissions  produced  in  a  three  species  aurora  are 
another  model  result  that  needs  to  be  compared  to  observations. 


5.  IONOSPHERIC  PLASMA  TURBULENCE  AND  PARTICLE 

ACCELERATION  ASSOCIATED  WITH  DISCRETE-ARC 

PRECIPITATION 

Auroral  hiss  is  a  wave  phenomena  commonly  observed  in  the  acceleration  region 
above  the  auroral  zone,  associated  with  inverted-V  electron  precipitation.  The  frequencies 
of  these  waves,  extending  from  the  lower-hybrid  resonance  up  to  the  plasma  frequency, 
fall  in  the  VLF  range,  and  their  dispersion  characteristics  have  1^  them  to  be  identified  w  ith 
waves  on  the  whistler  resonance-cone  plasma  dispersion  surface.  The  electric  field 
intensities  of  these  waves  are  occasionally  large,  and  this  turbulence  has  been  observed  to 
be  correlated  with  the  energy  flux  of  transversely  accelerated  ions  observed  in  the  topside 
ionosphere. 

A  mechanism  by  which  wave-particle  interaction  of  ions  with  this  intense 
turbulence  near  the  lower-hybrid  frequency  leads  to  ion  conic  fonnation  was  suggested  by 
Chang  and  Coppi;  the  turbulence  is  excited  by  the  instability  of  the  ambient  plasma 
combined  with  the  auroral  electrons  accelerated  to  drift  along  the  geomagnetic  field  by  the 
field-aligned  potential  drop  in  the  suprauroral  region.  Wave-particle  interaction  of  the 
ambient  ion  population  with  the  turbulence  near  the  lower-hybrid  frequency  leads  to  ion 
acceleration  nearly  perpendicular  to  the  field  line,  which  is  followed  by  the  adiabatic  folding 
of  velocities  as  the  ions  mirror  auid  travel  up  the  field  line,  to  create  the  conic  vekxity 
distribution. 

Plasma  simulations  by  us  have  shown  this  mechanism  to  be  effective  not  only  for 
the  transverse  acceleration  of  ions,  but  also  for  the  acceleration  of  the  ambient  electrons  in 
the  directions  pmallel  to  the  geomagnetic  field.  F.lcctron  vekxity  distributions  are  found  to 
have  enhanced  fluxes  of  energetic  electrons  in  both  directions  with  respect  to  the 
geomagnetic  field,  a  form  which  is  indeed  observed  in  conjunction  with  ion  conics,  where 
it  is  called  counterstrearning  electrons.  (We  refer  here  to  the  counterstreaming  electrons  of 
Tvpe  ’T”,  in  which  there  is  no  evidence  of  electron  acceleration  by  DC  pttrallcl  electric 
fields.) 
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5.1  Two  Dimensional  Plasma  Simulations 

In  the  interest  of  simplicity,  our  earlier  simulations  were  one  dimensuinal 
simulations  in  which  the  waves  could  propagate  in  only  one  direction  with  respect  to  the 
ambient  magnetic  field.  A  number  of  interesting  phenomena,  however,  can  be  realii'cd 
only  in  higher  dimensions,  where  a  range  of  propagation  angles  can  be  studied 
simultaneously.  For  example,  both  the  real  frequency  and  the  convective  linear  growth 
rates  of  waves  on  the  whistler  resonance-cone  dispersion  surface  are  a  function  of  the  angle 
of  propagation  with  respect  to  the  magnetic  field.  Mode  coupling  processes  cannot  be 
studied  in  complete  generality  in  one  dimension,  because  wave  vectors  there  arc 
constrained  to  be  aligned  in  the  same  direction.  Finally,  because  electrons  and  ions  interact 
with  different  efficiencies  to  waves  of  different  frequencies,  an  estimation  of  the  relative 
effectiveness  of  electron  and  ion  acceleration  requires  a  simulation  in  which  waves  of  an 
appropriate  range  of  frequencies  can  be  excited.  To  address  these  questions,  u  two 
dimensional  simulation  model  was  designed  and  a  series  of  simulation  runs  were 
performed. 

In  the  two-d  simulations,  strong  particle  acceleration  was  observed,  both  of 
electrons  and  of  ions.  Because  of  their  restricted  perpendicular  mobility,  electron 
acceleration  occurred  primarily  parallel  to  the  magnetic  field.  The  acceleration  caused 
energetic  tail  formation  on  the  electron  parallel-velocity  distribution,  both  in  the  direction  of 
the  beam  and  in  the  opposite  direction.  The  acceleration  in  the  direction  opposite  to  the 
beam,  as  well  as  the  initial  acceleration  in  the  beam  direction,  can  be  attributed  to  resonant 
interaction  with  short-wavelength  waves  excited  by  the  nonlinear  decay  of  the  long- 
wavelength  waves  excited  linearly  by  the  beam.  The  ions  tend  to  be  accelerated  in  the 
direction  perpendicular  to  the  magnetic  field  and  lighter  ions  are  accelerated  to  greater 
individual  energies  than  heavier  ions  are.  These  are  both  consequences  of  the  resonant 
nature  of  the  interaction  with  the  waves;  the  perpendicularly  propagating  waves  have  lower 
frequencies  and  phase  velocities  so  the  ions  can  more  easily  interact  with  them,  and  the 
lighter  ions  have  higher  thermal  velocities  and  will  be  more  likely  to  match  the  phase 
velocity  of  a  wave.  The  mass  dependence  of  the  velocity  diffusion  coefficient  also  favors 
the  acceleration  of  the  lighter  ion  species.  Given  the  observed  wave  spectrum,  the 
quasilinear  description  of  particle  acceleration  appears  to  be  adequate  to  account  for  both  the 
electron  and  ion  acceleration  observed. 

Much  of  the  interesting  physics  of  the  simulation  is  contained  in  the  evolution  of  the 
waves,  which  can  be  viewed  by  means  of  simulation  snapshots  of  electric-field  spectral 
densities  with  respect  to  frequency  and  wavevector.  In  the  classic  picture,  waves  are 
generated  at  relatively  long  wavelengths  by  the  linear  instability;  once  these  waves  become 
intense  enough,  they  succumb  to  a  variety  of  nonlinear  mode-coupling  processes  which 
attempt  to  thermalize  and  isotropize  the  spectrum.  This  process  can  be  documented  with 
the  snapshots  from  the  simulation.  The  mode -coupling  processes  also  aid  panicle 
acceleration  by  making  the  excited  waves  more  accessible  in  phase  velocity  to  the  low- 
energy  ions  and  electrons.  Spectral  analysis  of  the  temporal  evolution  of  the  potential  at 
one  position  helps  us  identify  the  waves  that  an  in-situ  observer  would  measure.  Doing 
this,  we  find  a  broadband  spectrum  of  frequencies,  ranging  from  the  lower  hybrid 
resonance  frequency  up  to  the  plasma  frequency.  Note  that  if  a  one-dimensional  simulation 
of  the  .situation  had  been  performed,  only  the  plasma  frequency  waves,  due  to  their 
stronger  linear  coupling  with  the  auroral  electrons,  would  have  been  produced.  The 
observations  of  auroral  hiss  indicate  that  the  spectrum  of  the  turbulence  is  in  fact  broadband 
like  the  results  of  the  two  dimensional  simulations.  Thus  we  see  the  importance  of  the 
generalization  to  a  simulation  model  with  more  degrees  of  freedom.  Because  ions  can 
resonantly  interact  only  with  low-frequency  wave.s,  their  demonstrated  acceleration  in  the 
two-d  simulations  and  in  space  reinforces  this  conclusion. 
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5.2  Comparisons  to  Observations 

Another  aspect  of  our  study  of  ion  acceleration  in  the  topside  ionosphere  involves 
the  analysis  of  observational  data  and  the  application  of  theoretical  imxlels  to  the  data  in 
order  to  test  and  refine  the  models,  'fhe  rocket  flight  MARIE  gave  us  what  is  probably  the 
best  set  of  data  illustrating  the  formation  of  ion  conics  through  the  wave-particle  interaction 
of  ions  with  intense  VLE  turbulence.  Ion  flux  detectors  recorded  energetic  ion  fluxes  near 
pitch  angles  of  90  degrees,  indicating  that  the  transverse  heating  of  the  ions  w  as  (xcurnng 
in  the  vicinity  of  the  rocket.  Mass  spectrometer  data  showed  that  the  energetic  ions  were 
predominantly  H-t-.  Simultaneously,  electric  field  measurements  indicated  an  intense 
auroral  hiss  spectrum,  with  its  low-frequency  cutoff  at  the  lower  hybrid  freijuency. 
Although  the  intensities  of  the  ion  flux  and  electric  field  correlated  strongly  with  one 
another  during  the  relevant  portions  of  the  flight,  it  was  still  necessary  to  establish  that  the 
w  aves  are  responsible  for  the  heating. 

To  do  this,  we  performed  a  series  of  Monte  Carlo  simulations  of  the  event,  using 
the  observed  electric  fields  to  heat  the  ambient  ion  population  to  determine  if  the  turbulence 
can  heat  the  ions  to  the  observ'ed  energies.  Pitch-angle  data  gave  us  some  information 
about  the  geometry  of  the  event  and  the  extent  of  the  heating  region,  which  was  less  than 
100  km  along  the  geomagnetic  field  line.  Without  information  about  the  wave-vector 
spectrum  of  the  waves,  we  had  to  rely  on  either  models  of  the  k-spectrum  or  simple  order- 
of-magnitude  estimates  of  it.  because  the  wave-panicle  resonance  is  a  velocity-dependent 
one.  For  these  calculations,  we  used  the  latter,  and  from  the  intensity  of  the  electric  field 
spectrum,  estimated  that  the  H+  heating  rate  was  about  0.8  eV/sec. 

From  the  Monte  Carlo  calculations,  we  found  that  the  turbulence  heats  the  H-*-  ions 
to  an  average  energy  of  about  7  eV,  starting  from  an  initial  temperature  of  1  eV,  with 
measurable  ion  flux  out  to  100  eV.  litis  is  smaller  than  the  observed  energy;  the  observed 
temperature  is  22  eV,  with  ion  flux  out  to  400  eV.  With  the  remaining  uncertainties  in  the 
model,  however,  the  fact  that  these  predicted  energies  are  comparable  to  the  ouscived 
energies  is  a  signal  of  success.  It  is  to  be  noted  especially  that  if  the  waves  of  the 
turbulence  are  being  absorbed  by  the  ions  as  they  are  being  excited  in  the  plasma,  the 
measured  level  of  the  electric  field  will  not  be  a  true  measure  of  the  strength  of  the  ion 
heating  rate. 


5.3  Future  Directions 

One  of  the  consequences  of  the  inhomogeneity  of  density  and  magnetic  field  in  the 
topside  ionosphere  is  expected  to  be  the  processes  by  which  the  characteristics  of  plasma 
waves  generated  in  one  region  of  space  change  as  the  waves  propagate  to  other  regions. 
To  examine  the  role  that  these  processes  play  in  determining  the  nature  of  the  pla.sma  waves 
and  particle  acceleration  observed  in  the  auroral  zone,  a  set  of  plasma  simulations  could  be 
designed.  By  launching  whi.stler  resonance-cone  waves  into  a  density  gradient,  one  could 
investigate  the  propagation  of  the  waves  into  the  region  where  their  frcqucucy  meets  the 
local  lower-hybrid  re.sonance  frequency  and  they  are  either  absorbed  directly  or  convened 
into  shorter  wavelength  waves  which  can  perhaps  be  more  efficiently  absorbed  by  the 
ambient  plasma  populations.  In  this  way,  it  is  expected  that  the  mode  conversion  process 
acts  as  an  aid  to  the  production  of  transverse  ion  acceleration  into  the  ion  conics  that  are 
observed  in  the  region. 
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MODELING  THE  IONOSPHERIC  E-REGION 


Over  the  years,  many  researchers  have  developed  E  region  models,  ou.'^sclvo 
included,  but  a  common  problem  has  been  a  tendency  of  the  models  to  underestiinaie  tfic 
electron  density  when  compared  to  observational  data.  While  a  variety  of  reasons  have 
been  considered,  we  have  found  in  the  literature  no  definitive  resolution  to  this  problem. 
As  recently  as  1988,  Rasmussen  et  al.  (1988)  discussed  this  problem  and  concluded  that 
little  improvement  was  possible  at  that  time  due  to  uncertainties  in  solar  fluxes,  reaction 
rates,  cross  sections,  etc.  More  recently,  Buonsanto  (1990)  compared  daytime  nuxlel 
results  with  ionosonde  data  over  a  full  solar  cycle,  a  variety  of  Millstone  Hill  radar  data, 
and  a  couple  of  ion  composition  models.  While  recognizing  the  same  problems  noted  by 
Rasmussen  et  al.  (1988),  Buonsanto  (1990)  had  reasonable  success  between  1  lO  and  )xn 
km  if  (1)  the  standard  EUV  solar  fluxes  were  increased  by  25-}i{)%  and  (2)  the  neutral 
densities  given  by  MSIS-86  were  decreased. 


6.1  Auger  Electron  Effects 

In  section  3.1,  we  described  the  role  of  solar  soft  X-rays  and  Auger  electrons  in 
contributing  to  the  high  altitude  energetic  photoelectron  flux.  In  this  section,  we  describe 
our  study  of  the  effects  of  solar  X-rays  and  Auger  electrons  on  the  Earth’s  E-region. 
Using  a  local  steady-state  model  of  the  Eanh’s  E-region,  we  have  found  that  the  importance 
of  the  Auger  electrons  depends  critically  on  the  intensity  of  the  solar  soft  X-rays  and  the 
solar  zenith  angle  (SZA).  At  large  SZAs  during  solar  maximum,  we  found  the  most 
dramatic  effects  from  the  Auger  electrons,  whereas  during  solar  minimum  and  at  small 
SZAs,  the  effects  are  much  reduced.  The  quantities  most  affected  are  those  that  depend 
directly  and  simply  on  the  photoelectron  flux,  such  as  airglow  volume  emission  rates  and 
electron  second^  production  rates.  On  the  other  hand,  quantities  with  less  sensitivity  to 
the  photoelectron  flux,  such  as  the  electron  and  ion  densities,  are  less  affected.  The  largest 
uncertainty  in  evaluating  the  role  of  the  Auger  electrons  comes  from  our  lack  of  knowledge 
concerning  the  intensity  of  the  solar  soft  X-rays,  a  region  of  the  solar  spectrum  that  is 
highly  variable  and  inadequately  measured.  In  an  attempt  to  evaluate  the  most  appropriate 
solar  fluxes  to  use,  we  surveyed  the  available  measured  solar  fluxes  and  the  scaling  models 
derived  from  them.  One  constraint  on  the  solar  flux  might  be  made  by  requiring  that  an  E- 
region  model  using  the  flux  produce  an  E-layer  density  which  is  in  agreement  with 
empirical  E-Iayer  models,  which  are  based  on  observations  of  E-region  densities.  With  this 
in  mind,  we  made  a  series  of  calculations  of  the  E-region  electron  density  profile  using  the 
Hinteregger  reference  flux  SC#21REFW  and  the  associated  FI 0.7  algorithms.  We  found 
that  the  E-region  peak  density  predicted  by  our  theoretical  E-region  model  is  consistently 
less  than  that  predicted  by  an  empirical  E-region  model.  To  get  agreement  in  the  solar 
minimum  case  (F10.7  =  77),  a  scaling  factor  of  3.4  was  applied  to  the  high  energy  solar 
flux  (E  <  88  eV).  While  not  small,  this  scaling  of  the  flux  is  reasonably  consistent  with  the 
findings  of  other  workers.  One  difficulty  in  this  approach  is  that  the  solar  flux  is  not  the 
only  input  of  theoretical  E-region  models  which  is  known  to  be  subject  to  large  uncertainty 
or  natural  variability.  In  particular,  consideration  of  the  role  of  NO  in  the  model  led  us  to 
develop  an  improved  characterization  of  the  NO  altitude  profile. 


6.2  Global  Modeling 

In  order  to  complement  our  ability  to  model  the  F  region  of  the  Earth's  ionosphere 
on  a  global  scale,  we  extended  what  was  an  E  region  daytime  model  to  include  the 
nighttime  E  region.  This  involved  implementing  a  nocturnal  source  of  ionization  similar  to 
that  in  our  F  region  model.  For  the  E  region,  we  included  the  production  of  N2-i-  and  O2+ 
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as  well  as  0+,  and  we  included  photons  at  a  tounh  wavelength  ( 1026).  In  order  to  handle 
the  large  solar  zenith  angles  around  sunset  and  sunrise,  our  calculation  ol  neutral  column 
density  was  completely  rewritten,  ft-eviously,  we  had  used  the  secant  approsuiiation  (or 
angles  less  than  80°  and  used  Chapman's  grazing  incidence  approximation  tor  angles  up  to 
90°.  While  the  Chapman  approximation  can  be  used  for  angles  above  90°,  our  approach 
had  the  problem  of  causing  a  discontinuity  when  we  switched  from  the  secant  to  the 
Chapman  approximation.  This  was  true  even  when  we  shifted  the  cross  over  point  to 
lower  angles.  Our  new  version  does  line  of  sight  numerical  integration  at  angles  above 
23.5°  and  uses  the  secant  approximation  below,  giving  a  highly  accurate  and  smoothly- 
behaving  column  density.  In  addition  to  a  nighttime  capability,  the  nnxlcl  can  now  be  run 
over  a  grid  of  longitudes,  latitudes,  and  times. 


6.3  Future  Directions 

As  mentioned,  in  our  model  we  have  scaled  the  solar  soft  X-ray  flux  above  the 
values  given  in  published  reference  spectra.  While  this  has  helped  bring  the  modeled  FoE 
in  agreement  with  observations,  between  the  E  layer  and  F2  layer  peaks,  we  still  see  the 
tendency  to  underestimate  the  observations.  Future  work  should  involve  a  study  of  the 
sensitivity  of  the  model  to  variations  in  neutral  density,  electron  temperature,  recombination 
rates,  and  0+(2D)  chemistry.  The  purpose  would  be  (1)  to  see  if  you  would  have  any 
success  with  the  scalings  suggested  by  Buonsanto  and  (2)  to  see  what  fine  tuning  of  the 
mtxiel  will  reproduce  reasonable  E  regions. 

To  make  funher  progress  on  the  role  of  the  solar  soft  X-rays,  one  could  do  a 
detailed  analysis  of  coincident  data  from  the  Low  Altitude  Plasma  Instrument  on  DE-2  and 
ground  based  digisondes.  In  such  a  study,  you  would  use  a  measured  electron  spectrum  to 
infer  a  .solar  flux  and  then  use  that  flux  to  calculate  an  E-region  electron  profile.  That 
profile  would,  in  turn,  be  compared  to  a  digisonde  measured  profile. 


7 .  MODELING  THE  MIDLATITL'DE  F-REGION 
7 . 1  Model  Modifications 

Over  the  life  of  the  contract,  we  have  made  a  series  of  modifications  to  our 
theoretical  0-t-  F-region  model.  The  first  modification  was  implementing  the  horizontal 
wind  model  (HWM87)  of  Fledin  el  al.  (1988)  into  our  code.  This  is  an  empirical  neutral 
wind  model  based  on  wind  data  obtained  from  the  Atmosphere  Explorer  E  and  Dynamics 
Explorer  2  satellites.  With  HWM87,  we  could  then  include  zonal  winds  in  our  mtxiel,  but 
to  take  advantage  of  this,  we  needed  to  include  the  magnetic  field  declination.  From  earlier 
work  (Rush  ei  al.,  1984),  we  knew  that  including  the  zonal  wind  and  magnetic  declination 
is  an  essential  element  in  modeling  the  longitude/UT  dependence  in  the  F  region 
ionosphere.  At  the  time  of  the  earlier  studies,  reliable  neutral  wind  models  for  use  in 
ionospheric  modeling  were  not  available.  Bui  with  models  available  such  as  HWM87  and 
the  Vector  Spherical  Harmonic  model  of  Killeen  et  al  (1987),  the  time  was  right  for 
studying  the  importance  of  the  ionospheric  effects  of  the  zonal  wind  and  magnetic 
declination. 

Two  other  inputs  in  our  E-region  mtxiel  are  the  Te  and  Ti  profiles  and  they  can  have 
a  strong  impact  on  a  model  s  ability  to  agree  with  measurements.  The  original  nuxlel  used 
an  e.mpirical  Te  and  Ti  mtxiel  based  on  Strobel  and  McElroy  (1970),  We  have  now 
upgraded  both  the  Te  and  Ti  mtxlels.  For  Te,  we  decided  not  to  u.se  the  model  as  found  in 
the  International  Reference  Ionosphere  (IRI)  but  to  implement  our  own  version  of  the 
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global  empirical  Te  model  of  Brace  and  'I'heis  (198 1 ).  This  nuKiel  gives  Te  at  3(X)  and  4(/{ ) 
km  as  a  function  of  season,  dip  latitude,  and  hKal  time.  We  constructed  an  entire  prot'iic 
through  a  combination  of  linear  interpolation  and  a  profile  shape  function  b;r  cd  <tn 
Millstone  Hill  radar  data  (Strobel  and  McElroy,  1970).  f  'or  the  ion  temperature.  \u*  turned 
to  ion  temperature  measurements  made  at  Millstone  Hill  fur  the  year  1964  (Hvans,  1467 1, 
From  ihis  data,  we  constructed  a  tabulation  of  Ti  at  4{K)  and  6(X)  km  as  a  function  of  season 
and  local  lime.  The  entire  profile  is  again  constructed  through  a  mixture  of  hnear 
interpolation  and  a  profile  shape  function. 

The  next  changes  involved  implementing  new  diffusion  and  chemistrv'  rates  so  as  to 
be  consistent  with  existing  low  and  high  latitude  models.  Finally,  we  improved  our 
daytime  .source  of  0+  in  our  F  region  model.  Originally,  the  source  was  calculated  using  a 
very  crude  approximation.  The  solar  flux  was  characterized  by  a  single  number  and  the 
absorption  of  photons  in  the  atmosphere  was  calculated  using  a  single  cross  section. 
Funher.  there  was  no  explicit  source  due  to  secondary  ionization,  but  rather  it  was  ptnuiv 
approximated  by  not  allowing  0+(2D)  to  be  convened  to  N2+.  The  result  was  that  while 
computationally  fast,  these  approximations  .severly  underesiimaicd  0+  production  at  lower 
altitudes.  When  calculating  the  full  electron  density  profile,  we  could  compensate  by  using 
at  lower  altitudes  the  0+  density  from  our  E  region  model.  This  meant  that  in  order  to  get  a 
full  0+  profile,  we  had  to  run  both  models.  On  the  other  hand,  especially  for  global 
modeling,  we  wanted  to  be  able  to  use  the  F  region  model’s  results  directly  without 
depending  on  the  E  region  model.  Unfonunately,  the  complete  calculation  of  the  0+ 
source,  which  we  u.se  in  our  E  region  model,  is  too  slow  for  u.se  in  the  F  region  model. 
Thus,  we  have  developed  an  0+  source  routine  which  is  sufficiently  accurate  but  fast 
enough  for  our  purposes.  It  divides  the  solar  flux  and  various  cross  sections  into  i ! 
energy  bins  and  includes  a  calculation  of  the  photoionization  of  O,  a  calculation  of  the 
dissociative  photoionization  of  O2,  and  an  estimate  of  the  electron  impact  ionization  of  O 
based  on  the  work  of  Richards  and  Torr  (1988).  This  gives  a  much  improved  bottomside 
0+  profile  and  allows  the  calculation  of  the  total  electron  density  profile  by  simply  adding 
0+  and  H+  profiles  from  the  F  region  model  with  O2+  and  N0+  from  the  E  region  mcxlcl. 


7,2  Theoretical  Studies 

In  this  section,  I  describe  two  studies;  (1 )  the  role  of  the  zonal  wind  and  magnetic 
declination  and  (2)  a  sensitivity  study.  I  begin  with  the  zonal  wind  study.  The  importance 
of  including  the  effects  of  neutral  winds  in  any  solution  of  the  continuity  etjuadon  is  that  in 
the  F  region  the  neutral  wind  moves  ionization  along  geomagnetic  field  lines.  When 
considering  this  interaction,  it  is  usually  the  meridional  component  of  the  neutral  wind  that 
gets  all  the  attention.  This  is  bee,  use  the  magnetic  field  is  roughly  aligned  w  ith  the 
meridional  wind  such  that  the  equatorward  meridional  wind  moves  ionization  up  in  altitude, 
and  the  poleward  meridional  wind  moves  ionization  downward.  However,  this  alignment 
is  not  exact.  Rather,  the  magnetic  field  does  have  a  declination  and  thus,  the  zonal 
component  of  the  neutral  wind  can  have  a  projection  along  the  field.  The  result  is  that  the 
zonal  wind  imparts  an  upward  or  downward  drift  to  the  ionization  depending  upon  the 
relative  direction  between  the  v'ind  velocity  and  the  declination  of  the  geomagnetic  field. 
The  interaction  between  the  zonal  wind  and  the  geomagnetic  field  gives  rise  to  very 
different  diurnal  variations  of  the  0+  density  profile  at  eastward  and  westward  declinations. 
It  is  expected  that  the  clearest  signature  of  this  effect  can  be  found  by  looking  at  the 
longitude  dependence  of  the  diurnal  behavior  of  the  0+  profile.  While  this  effec  t  was 
studied  by  Rush  et  al.  (1984),  they  used  a  fairly  simple  neutral  wind  as  well  as  a  neutral 
atmosphere  with  no  longitudinal  dependence.  By  using  the  horizontal  wind  model 
(HWM87)  of  Hedin  et  al.  (1988),  along  with  the  MSIS86  neutral  atmosphere,  we  can  do  a 
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more  through  study  of  diurnal  variations  at  different  longitudes  and  hence,  a  more  cotnpieie 
assessment  of  the  importance  of  the  zonal  wind  and  magnetic  declination. 

For  our  first  sets  of  calculations,  we  used  a  lilted  dipole  magnetic  field  mixiei  to 
specify  the  magnetic  inclination  and  declination.  We  made  some  initial  calcuiaiions  at  two 
locations  that  had  different  longitudes  (12  east  and  210  east)  but  had  the  .same  geographic 
(40  nonh)  and  geomagnetic  (40.8  nonh)  latitudes.  In  this  way,  we  could  examine  the 
differences  that  arise  between  kxtations  with  the  same  solar  forcing  and  inchnations  bui 
different  longitudes.  We  found  that  there  were  differences  in  the  0+  densiiics  ai  these  two 
locations.  For  example,  an  examination  of  the  diurnal  variation  of  the  total  ion  content 
(TIC)  found  that  at  1(K)0  ItK'al  time  (LT)  the  TIC  at  210  east  is  50  percent  higher  than  the 
KXX)  LT  TIC  at  12  east.  At  night,  this  relationship  was  seen  to  reverse  and  the  22vX)  LT 
TIC  at  12  east  was  three  times  larger  than  that  at  210  east. 

Our  next  step  was  to  expand  our  calculations  to  a  more  global  scale.  The  first 
global  calculations  were  for  the  northern  hemisphere  at  two  geomagnetic  longitudes  (90.98 
(E)  and  269°  (E))  and  over  a  range  of  geomagnetic  latitudes  (70(N)  to  26(N)).  The  tilted 
dipole  model  for  the  Earth's  magnetic  field  gave  declinations  at  these  two  longitudes 
opposite  in  sign  and  ranging  in  magnitude  from  37.4  to  13  degrees.  At  a  given 
geomagnetic  latitude,  each  longitude  had  the  same  inclination  and  geographic  latitude. 

The  calculations  done  for  just  two  locations  had  given  us  our  fir.st  examples  of 
longitudinal  differences,  but  as  we  began  the  global  calculations,  we  came  to  appreciate  that 
a  longitudinal  dependence  does  not  simply  reflect  differences  in  the  magnetic  declination. 
Rather,  within  the  context  of  the  mcxiel  at  a  fixed  latitude,  the  longitudinal  differences  in  the 
diurnal  behavior  are  due  to  differences  in  declination,  longitudinal  dependencies  in 
FIV»'M87,  and  longitudinal  dependencies  in  the  MSLS86  neutral  atmosphere  model.  To 
separate  these  three  effects,  we  made  three  sets  of  calculations.  The  first  set  consisted  of 
full  calculations  of  the  diurnal  behavior  of  the  0+  density  along  the  two  magnetic 
longitudes  given  above.  'Fhe  second  set  were  the  same  calculations  e.xccpt  the  declinations 
were  set  to  zero,  which  effectively  turned  off  the  zonal  wind.  In  the  third  set,  we  repeated 
the  calculations  except  both  components  of  the  neutral  wind  were  set  to  zero. 

From  the  first  set,  we  find  variations  with  longitude  at  all  latitudes  similar  to  what 
we  saw  at  one  latitude.  The  general  features  are  the  same  from  latitude  to  latitude  though 
the  quantitative  details  do  vary  markedly.  These  differences  in  detail  rellect  the  latitudinal 
dependence  of  the  declination  and  inclination  as  well  as  the  latitudinal  dependencies 
contained  in  the  HWM87  wind  model  and  the  MSIS86  neutral  atmosphere  model. 

When  we  focus  on  a  given  latitude  and  compare  the  results  from  ihe  different  scin. 
the  different  longitudinal  effects  become  apparent.  Comparison  of  set  1  and  set  2  illustrates 
the  role  of  the  zonal  wind  and  magnetic  field  declination.  We  find  that  at  30  degrees 
latitude  from  0800  to  13(X),  most  of  the  variation  with  longitude  arises  from  differences  in 
declination  and  any  longitudinal  dependencies  in  the  HWM87  zonal  winds.  From  13(X)  to 
Of/X"),  the  variation  with  longitude  is  panially  effected  by  zonal  winds  and  declination  but  a 
large  percentage  of  the  variation  is  due  to  longitudinal  variations  of  the  meridional  wind 
and/or  longitudinal  variations  in  the  MSIS86  neutral  atmosphere. 

Comparisons  of  set  1  and  set  3  show  that  the  neutral  atmosphere  plays  no  role  in 
the  longitudinal  variations  at  night  and  for  this  case,  at  30  degrees,  a  small  role  around 
noon.  Set  3,  with  no  neutral  winds,  does  illustrate  the  general  point  that  the  neutral  winds 
do  make  a  dramatic  difference  in  the  calculated  0+  densities.  A  direct  examination  of  the 
neutral  winds  at  the  two  longitudes  shows  that  they  are  consistent  with  the  variations  with 
longitude  seen  in  the  0+  densities.  Again,  these  features  are  seen  at  all  latiiudcs  to  differing 
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degrees.  From  these  initial  results,  we  see  that  the  diurnal  variation  of  the  F2  regain  is 
strongly  effected  by  the  zonal  wind  and  magnetic  declination,  as  well  as  the  longituihnai 
and  latitudinal  variations  in  the  meridional  wind. 

We  expanded  on  this  initial  work  by  including  a  more  thorough  parameter  study. 
That  is,  in  the  initial  calculations,  we  had  considered  just  one  set  of  geophysical  conditions. 
Now,  we  made  calculations  foi  three  seasons  (March  equinox,  June  solstice,  and 
December  solstice)  and  three  levels  of  magnetic  activity  (Ap  =  5.  11,  32).  We  also  did  sets 
of  five  calculations  for  each  combination  of  geophysical  conditions.  In  the  first  calculation 
of  a  set,  we  calculated  the  diurnal  behavior  of  the  0+  density  with  the  longitude  dependence 
in  both  HWM87  and  MSIS86  turned  off  and  the  zonal  wind  set  to  zero.  In  the  next,  we 
had  the  longitude  dependence  in  MSIS86  turned  on.  The  third  calculation  had  the  longitude 
dependence  in  MSIS86  turned  on  and  included  a  longitudinally  independent  zonal  wind. 
The  founh  had  the  zonal  wind  set  back  to  zero  but  included  a  longitude  dependent 
meridional  wind  from  HWM87  along  with  the  longitudinally  dependent  MSIS86.  Finally, 
the  fifth  calculation  had  all  the  longitude  dependencies  turned  on  and  included  the  zonal 
wind.  In  total,  90  simulations  of  the  northern  midlatitudes  were  conducted. 

What  we  found  was  consistent  with  our  earlier  results.  All  three  factors  described 
above  played  a  role.  In  general,  the  critical  sources  of  longitude  differences  were 
zonal/declination  differences  and  the  longitude  dependence  of  the  meridional  wind  given  by 
HWM87.  The  observed  differences  in  the  daytime  morning  were  dominated  by  zonal/ 
declination  effects,  while  the  nighttime  differences  were  driven  by  a  mix  of  zonal/ 
declination  effects  and  longitude  dependent  meridional  winds.  While  the  details  varied,  the 
above  description  was  true  at  all  latitudes,  seasons,  and  magnetic  activities  that  we 
considered.  The  effects  tended  to  be  larger  as  you  approached  the  equator.  They  were  also 
larger  for  more  active  magnetic  conditions.  The  largest  effects  were  in  the  summer  and  the 
smallest  in  winter.  In  all  cases,  in  the  daytime  morning,  the  91°  (E)  longitude  calculation 
produced  smaller  densities  than  the  269°  (E)  calculation,  and  at  night  the  reverse  was  true. 

In  our  second  study,  we  made  a  series  of  calculations  where  we  individually  varied 
many  of  the  inputs  of  the  model.  Two  of  these  calculations  involved  variation  of 
microscopic  parameters:  the  0-0+  collision  cross  section  and  the  fraction  of  0+(2D)  that 
converts  to  N2+.  The  remaining  calculations  involved  variations  of  the  different 
geophysical  inputs:  the  electron  temperature,  the  atomic  oxygen  density,  the  neutral 
molecular  densities,  the  solar  flux,  and  the  neutral  wind.  All  of  these  parameters  effect  the 
F2  peak  density  but  only  the  molecular  densities  and  the  neutral  wind  have  any  significant 
effect  on  the  F2  peak  altitude. 

There  has  been  some  debate  in  the  recent  literature  as  to  whether  the  accepted  value 
for  the  0-0+  collision  cross  section  should  be  increased.  Increasing  the  cross  section  by 
the  suggested  factor  of  1.7  leads  to  just  a  5%  change  in  the  F2  peak  density.  Going  from 
one  extreme  of  having  all  the  0+(2D)  conven  to  N2+  to  having  none  of  it  convert  produces 
a  35%  increase  in  the  F2  peak  density.  Turning  to  the  electron  temperature,  we  find  that  a 
30%  scaling  leads  to  just  a  5%  effect  on  the  peak  density.  For  both  the  atomic  oxygen  and 
the  entire  solar  flux,  the  modeled  F2  peak  density  shows  a  basically  linear  dependence. 

As  mentioned,  the  molecular  neutral  densities  effect  both  the  F2  peak  density  and 
altitude.  A  30%  variation  in  N2  and  €>2  leads  to  a  15  to  20%  variation  in  the  peak  density 
and  a  10  km  variation  in  the  peak  altitude.  But  of  greatest  interest  is  the  neutral  wind 
because  it  not  only  effects  both  the  peak  height  and  density  but  is  probably  the  geophysical 
input  that  is  most  variable  and  least  known.  For  example,  a  50%  reduction  of  the  neutral 
wind  from  the  Hedin  et  al.  (1988)  wind  model  (HWM87)  can  lead  to  a  12%  increase  in  the 
peak  density  and  a  15  km  increase  in  the  peak  height.  Further,  if  we  use  the  vector 
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spherical  harmonic  (VSH)  wind  model  of  Killeen  ei  al.  (19X7)  rather  ui.ui  hVvMh/,  we  eet 
a  peak  density  22%  larger  than  that  from  HWM87  and  a  peak  altitude  km  higtict,  Au 
examination  of  the  daytime  meridional  winds  from  these  two  models  snows  significant 
differences  which  lead  to  the  differences  seen  in  the  F2  peak  paramciers.  In  principle, 
these  two  models  are  describing  the  neutral  wind  field  for  the  same  geoptiysicai  conditions. 
These  differences  in  results  illustrate  two  points;  t ! )  an  ionoispheric  rnixie!  is  very'  sensitive 
to  uncertainties  in  the  neutral  wind  and  (2)  present  w-ind  models  can  predict  sigmiicantly 
different  winds  under  identical  geophysical  conditions.  Our  conclusion  is  that,  of  all  the 
inputs  to  our  ionospheric  model,  the  neutral  wind  is  the  most  likely  candidate  as  a  major 
source  of  error. 


7.3  Comparisons  with  Observations 

As  part  of  our  development  of  a  daytime  E  and  E  region  ionospheric  model,  we 
analyzed  a  series  of  experiments  involving  obserx-ations  by  the  Millstone  1  lill  Radar  facility 
coincident  with  daytime  overflights  by  a  DMSS  satellite.  Observations  were  made  from 
February  to  October  1989  involving  the  F9  satellite  {0930  Itx'al  time)  and  the  F8  .satellite 
(1800  local  time).  The  early  experiments  consisted  of  approximately  four  hours  of  radar 
elevation  scans.  Tliis  produced  a  time  history'  of  the  electron  density  profile  (EiDP)  before, 
during,  and  after  a  satellite  overilight.  In  later  experiments,  by  considering  only  very  close 
overflights,  the  Millstone  zenith  antenna  was  used  in  the  high  resolution  OA.S1S  {Optimal 
Analysis  of  Signals  from  Incoherent  Scatter)  mode  and  the  steerable  antenna  was  used  to 
obtain  vector  drift  measurements. 

We  performed  an  analysis  on  ten  experimettts  using  our  F  and  E  region  models  and 
our  generalized-inverse  least  squares  (GILS)  algorithm.  The  GILS  algorithm  w'as 
developed  in  earlier  studies  on  the  use  of  daytime  airglow  for  remotely  sensing  the  EDP. 
In  this  analysis,  we  ran  our  ionospheric  models  in  both  their  ab  initio  and  one  con.straint 
mode.  In  the  one  constraint  mode,  we  use  the  measured  ion  density  at  the  satellite  altitude 
as  an  indirect  constraint  on  the  geophysical  input  parameters  of  our  models.  By  using  the 
GILS  algorithm  to  find  the  most  likely  set  of  input  parameters  for  which  our  model  matches 
the  observational  constraint,  we  calculated  the  most  likely  electron  density  profile  (EDP), 
This  one  constraint  result  can  be  compared  with  the  ab  initio  result  and  the  radar  profile,  as 
well  as  the  EDPs  predicted  by  the  International  Reference  Ionosphere  (IRI)  and  the 
Ionospheric  Conductivity  and  Electron  Den.sity  (ICED)  model.  In  general,  we  found  that 
our  models  in  either  of  their  mixles  perform  better  than  IRI  and  are  comparable  to  ICED. 

While  we  had  made  a  few  comparisons  with  Millstone  radar  data  and  DMSP  in-situ 
data,  most  of  our  time  had  been  spent  on  modifying  our  model  and  examining  the 
theoretical  results  that  it  produced.  Thus,  we  saw  a  real  need  to  fcx-iis  on  whether  our 
results  were  realistic  and  if  our  results  agreed  with  other  models.  Our  approacii  was  to  sec 
what  "tuning"  of  our  midlatitude  F  region  model  was  necessary  in  order  to  satisfy  two 
criiena;  1 )  The  midlatilude  model  at  its  low  latitude  boundary  should  agree  with  a  low' 
latitude  model  at  its  high  latitude  boundary  and  2)  The  model  should  produce  realistic 
densities  as  compared  to  monthly  mean  data  or  empirical  models  of  monthly  means.  Our 
interest  was  not  in  attempting  to  model  individual  days  but  rather  to  see  if  using  "typical"  or 
average  inputs  in  the  model  would  produce  a  "typical"  or  average  ionosphere. 

Our  first  attempt  to  deal  wdth  both  of  these  issues  was  a  calculation  of  the  diurnal 
behavior  of  the  F2  peak  at  solar  maximum  (!'i().7  =  180)  in  March  for  Rome,  Italy. 
Initially,  we  had  differences  between  our  niidlatitude  model  and  a  low  latitude  model, 
'fhese  differences  were  found  to  be  caused  by  various  small  "hugs"  in  the  Uw'.  latitude 
code.  When  they  were  corrected,  the  two  models  agreed  to  within  a  few  percent. 
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Unfortunately,  while  we  achieved  agreeineni  between  intKiels.  agreenicru  wii'i  ror.  ia.mie 
data  was  problematical.  While  data  from  March  13-18  gave  midday  12  peak  dcii  aucN  of 
2. 1-2.3  X  10^  ciir-\  our  calculations  initially  ran  around  1.2  x  10'*  cm  •  By  niakiiic 
adjustments  in  the  neutral  winds  and  the  O-O-r  diffusion  coefficient,  we  were  anie  to 
increase  our  midday  densities  up  to  1,5  x  10^  cm'^.  But  obviously,  the  model  siiii 
underestimated  the  data.  We  were  also  awtu'e  that  other  modelers  had  similar  piof!k:iv.s 
modeling  for  soktr  maximum  conditions.  With  these  problems  in  mind,  we  decided  for  rhe 
time  being  to  focus  on  solar  minimum  conditions. 

For  the  solar  minimum  case,  we  focused  on  coinpari^'>ns  betvcs.,:!,  modci 
calculations  and  monthly  averaged  F2  peak  density  measureniems  made  at  Wallo'is  Isl.eiu, 
Virginia  in  March,  June,  September,  and  December  1986.  Our  initial  calculations  produced 
daytime  densities  that  were  in  qualitative  agreement  with  the  data,  but  for  three  of  the 
months  (March,  September,  and  December),  our  nighttime  densities  were  much  lower  than 
w'hat  was  measured.  It  was  only  in  our  June  calculation  that  rea.sonable  nighttiinc  densities 
were  obtained.  The  largest  differences  were  found  in  the  December  case  w  ith  calculated 
0-t-  densities  at  night  of  around  102  cm‘^  compared  with  measured  densities  of  around 
cm'^. 


Another  problem  with  the  December  case  was  a  numerical  instability  at  the  higher 
altitudes  (600  -  700  km)  that  began  a  bit  before  midnight  and  grew  until  photoionization  at 
sunrise  began  to  build  up  the  0+  densities.  At  first,  this  seemed  to  happen  only  when 
using  the  Brace  and  Theis  electron  temperature  model,  but  subsequent  calculations  showed 
that  the  instability  could  also  appear  when  using  our  older  temperature  model.  We  came  to 
believe  that  the  problem  was  related  to  the  extremely  small  densities  that  we  were  praiucing 
at  solar  minimum  nighttime.  A  check  of  earlier  calculations  uncovered  other  cases  w  here 
smaller  instabilities  could  be  seen  in  the  calculations.  These  cases  also  involved  small 
nighttime  densities. 

In  examining  the  December  case,  we  noted  that  one  reason  for  the  small  densities  at 
night  was  that  the  meridional  wind  did  not  turn  equatorward  until  well  after  sunset.  Using 
the  other  available  wind  model,  the  VSH  model  of  Killeen  et  al.  (1987),  gave  different 
resulLS  in  detail  from  those  gotten  using  HWM87  but  still  did  not  maintain  the  nighttime 
densities.  While  in  individual  cases,  you  could  argue  that  the  actual  neutral  wind  can  be 
quite  different  from  the  neutral  wind  models  we  still  had  to  deal  with  those  eases  where  the 
model  winds  v/ere  reasonable.  To  help  maintain  the  nighttime  density,  we  decided  to  add  a 
nocturnal  source  of  ionization  to  the  model.  We  used  an  approach  following  Knudsen  et 
al.  (1977)  of  simulating  the  ionization  due  to  scattered  EUV  photons.  We  assume  that  the 
scattered  EUV  photons  consist  of  three  lines  at  834,  584,  and  3()4A.  We  then  calculate  the 
photoionization  rates  as  if  these  three  lines  are  incident  at  the  top  of  the  ionosphere  with 
zero  zenith  angle  with  some  assumed  intensity.  Using  intensities  based  on  Knudsen  et  al. 
(1977)  and  Chakrabarti  et  al.  (1984),  we  found  that  the  nocturnal  source  helps  to  keep  the 
nighttime  den.sities  from  going  extremely  low  but  has  little  effect  much  above  a  density  of  a 
few  times  10^  cnr^. 

An  examination  of  the  literature  showed  that  maintaining  the  winter  nighttime 
ionosphere  has  been  a  problem  for  modelers  for  a  number  of  years.  While  there  has  been 
much  debate  as  to  what  is  the  necessary  physics  to  explain  the  observations,  it  seemed  to  us 
that  recent  modeling  indicated  that  including  hydrogen  ions  and  the  coupling  between  the 
ionosphere  and  the  plasma.sphere  will  be  necessary  for  any  successful  mtxieiing  of  the 
winter  nighttime.  Our  model  did  not  have  this  capability;  thus  we  decided  to  put  in  an  ad- 
hoc  source  to  maintain  our  nighttime  densities.  From  our  initial  experiment.s  w  ith  this  ad 
hoc  source,  we  found  that  when  the  nighttime  densities  are  maintained,  the  numerical 
instability  problems  do  not  appear. 
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VVe  next  turned  ti'  ns.jXHnuni  s.x>iidaii>ii>  by  m- -xieliiy  due  -oude  Jut.i  ai 
Waiktp',  Isiand  in  I-ebruary  ^>1  IvSd.  l-,;r  the  nioiuti  o:  heiiruarv,  \se  in  itoin:.-,;  a  day 
day  sitnulation  of  the  cMtue  na'ini'.  In  eompanny  taui'  results  vsuh  th.  VS  .i'din  -  Islaiui  data. 
vvC  examined  jUst  the  1-2  peak  densuy  and  hesent.  When  ue  etnnpaied  o.n  avciaye  mosicl 
results  to  the  average  davtintc  Ix-hasiuroi’  the  !'2  pe.'k  densuy.  ue  touad  ih.u  u,e  nuH^ei  did 
a  fair  job  in  that  it  was  pretty  imieh  within  a  standard  deviation  oS  tf;e  uio-cr\ed  na  aii  b!i! 
on  the  other  hand,  the  (laymste  p.>r!!t.!ji  uh  the  valeuiated  duirii.;i  sha;  e  ut  tfa  12  ne.tk 
density  svas  consistently  ditlesa-iu  faun  what  v-as  ofi-ienved.  I  he  I'i*  .ervaUons  i on  asteniiy 
had  the  daytime  maxiniutii  at  around  10. dO  ahcretis  the  mode!  iiad  tiic  tnasmuim  arotnui 
14:30.  Tills  was  true  on  a  daily  basis  as  weii  as  wuii  die  monthiv  .octaye.  .\  uuantitati vc 
comparison  of  the  averages  showed  that  at  i<..t:3l,‘.  tlie  mts  -iel  ‘Wits  ....'ll  i  jv»vs .  at  notin  it  was 
!  ITf  low,  and  at  14:. M)  it  was  within  IfT  of  me  data.  Wc  siutulu  note  if.ai  we  did  n<.>t  sec 
the  factors  of  two  differci...:es  that  have  been  seen  in  other  cienp-arisi'ns. 

When  we  compared  the  modeled  and  measured  i-2  peak  hcighi.s.  similar  results 
svere  seen.  That  is,  the  average  'ucxiel  results  were  within  a  sian.iaui  devuiuou  of  the 
monthly  averages!  data.  But  a  closer  examination  ■  "‘'e  entire  month  siiowed  that  from 

10:30  to  16:30,  the  nuKlcI  was  consistently  do'*'”  n'  •  sn  me  data,  predicting  an  1-2  i>e.ik 
heiglii  10  to  20  km  smaller  than  observ 

'I'he  next  series  of  si  ar  maximum  calculations  insob.csl  '^\cn  stations  and  four 
seasons.  The  stations  were  ly.mnui  '  islsiid.  and  MiUsti-.nc  Ihii  m  the  .Amencan 

longitude  sector;  .Slough.  .Moscow,  ana  f.,:iiiiingrau  from  the  iOirripean  sector;  .md 
Xinxiang  from  the  Asian  senior.  These  were  designed  as  monthly  mean  calculations  by 
using  the  monthly  averaged  h  lO.?  attd  .Ap.  The  l  ifteenih  of  each  month  vs  as  used  in  a!! 
eases  and  the  months  eho.scn  were  12/00,  3/00.  6/9(i  and  10/K9.  wiiii  the  exception  of 
.Moscow ,  where  the  months  mcKieled  were  1/90.  ’Si'-K).  7/H9  and  I 

Unfortunately,  data  became  available  only  for  .Millstone  Hiii.  V.killops  i.shmd.  and 
Bennuda.  /\l  Millstone,  comparisons  have  iven  made  for  Ociolver.  19,S9:  March.  1990; 
June,  1990;  and  December.  1990.  f-'or  the  December  case,  we  found  that  the  nnxleled  fulT 
at  nighttime  falls  below  the  scatter  in  the  data  at  1  to  2  .Mlu  below  an  .ivcrage  of  alx'ut  5 
■Mhz.  From  sunrise  until  KXK)  hours,  the  mrxJci  and  data  increase  together,  but  from  KKKI 
to  around  1200,  the  model  runs  beneath  the  monthly  scatter  again  at  about  1  to  2  Mh/, 
below  daytime  maximums  between  12  and  I  f  Mh/..  In  the  afternoon,  the  imxlel  falls  in 
the  lower  part  of  the  monthly  scatter  but  mvn  itig  into  evening,  u  falls  oil  a  bit  faster  than 
the  data.  Turning  to  the  March  case,  little  can  be  said  other  than  th.it  tlic  mnciel  tails  well 
w  ithin  the  enormous  monthly  scatter  of  the  data.  The  October  ease  is  siniihu  to  December 
except  that  the  monthly  sc;itier  is  large  enough  th.it  the  model  remains  wiiiiin  a.  Finally,  in 
June  the  observations  show  a  very  fiat  diurnal  curve  in  marked  contr.ist  te)  ilie  diurn.i! 
behavior  seen  at  other  seasons.  The  model  retJccts  this  change  m  darn.tl  shape  but  at  all 
limes  dearly  falls  at  the  high  extreme:-;  of  the  monthly  scatter  and  tluis,  consisiently  higher 
than  the  monthly  mean.  The  difference  is  in  the  1  to  2  Mh/  ’■aage  when  using  the  cJder 
0-+-/0  diffusion  rate  but  can  be  as  large  as  3. .3  .Mh/  when  using  tlv'  Burnside  Uk  tor  m  ihe 
ditfusion  rate, 


'I'he  companson.s  Ixtwccn  the  nuxlelcd  and  observed  hml2  give  a  similar  mixture  ol 
results.  In  all  four  cases,  the  model  falls  within  the  moninly  scatter  of  the  data.  For 
December,  the  agreement  with  the  mean  is  not  bad  in  the  daytime  but  .it  night,  the  model 
runs  about  50  km  higher  tlian  the  data.  In  June  we  sec  the  opposite,  the  nnxiel  is  around 
.50  krn  higher  from  0700  to  midnight  and  gives  fair  agrccmcnl  for  the  rest  of  the  day. 
.March  is  not  bad  o’ver  the  entire  tiay  but  again  the  scatter  is  large.  In  October,  there  is  a  bur 


scatter  but  the  nuxiel  does  trend  soineuhat  higher  than  ilu-  nican.  siuniai  icain  Ui  the 
December  case. 

For  the  other  locations  (Wallops  and  Bermuda)  similar  commeius  can  Ix’  made  In 
general,  we  do  best  in  the  winter  and  have  our  biggest  problems  in  the  summer. 

For  several  of  the  stations,  sve  were  surprised  by  the  results  from  around  June 
solstice  (May  -  July).  For  example,  the  calculation  done  Ujr  .Moscow  in  .May  sliowed  a 
diurnal  pattern  for  fof2  that  was  the  oppo.site  of  the  "texttxxik  "  pattern.  It  peaked  just  after 
2000  hours  local  tinte  at  around  1.^  Mhz,  stayed  flat  unul  02{)0  hours,  began  k,' decrease  to 
a  minimum  of  9  Mhz  just  after  1(X)0  hours,  stayed  flat,  and  finally  around  18(K)  hours, 
began  to  increase  towards  its  maximum.  hmf2  showed  a  large  change  going  trom  3(K)  km 
around  noon  to  over  480  km  at  around  22(X)  hours.  Similar  behavior  was  seen  in  Julv  but 
not  in  January  or  November.  These  results  are  in  disagreement  with  the  predictions  of  the 
International  Reference  Ionosphere  (IRI)  and  the  Field-Line  Interhemispheric  Plasma 
(FLIP)  model.  Our  initial  suspicion  that  this  result  wa»  due  directly  to  the  neutral  wind 
turned  out  to  be  correct.  Calculations  at  Bermuda  as  well  as  calculations  for  Moscow 
showed  that  .setting  the  wind  to  zero  significantly  modifies  both  the  fof2  and  hmil  giving 
diurnal  curves  that  are  more  "normal”.  An  examination  of  the  neutral  wind  being  used 
from  HWM87  showed  that  a  strong  meridional  wind  that  hardly  ever  turns  poleward  was 
responsible  for  our  surprising  results.  The  wind  used  in  the  FLIP  model  for  .Moscow  was 
deduced  from  the  hmf2s  predicted  by  the  IRI.  ITiis  wind  was  seen  to  be  radicady  different 
from  that  of  HWM87  and  accounts  for  the  differences  between  our  results  and  FLIP'S 
The  question  as  to  what  is  actually  going  on  at  Moscow-  is  unanswered  since  we  had  no 
Moscow  data  to  examine.  But  our  Bermuda  results,  while  high  compared  to  the  data,  did 
show  a  diurnal  shape  that  compared  well  with  data.  A  brief  examination  of  the  day-to-day 
Bermuda  data  showed  that  there  can  be  days  where  the  nighttime  densities  can  be  a  factoV 
of  2  larger  than  in  the  daytime. 


7.4  Future  Directions 

Present  midlatitude  ionospheric  models  have  been  successful  at  reproducing 
ionospheric  climate  (Sojka,  1989),  but  developing  models  that  can  reproduce  an  actual 
ionosphere  for  any  given  instant  will  require  much  more  extensive  comparisons  between 
models  and  observations.  Today,  we  are  in  a  position  to  make  progress  in  modeling  and 
testing  the  midlatitude  models  on  a  global  scale.  The  two  primary  reasons  that  this  is  true 
are  as  follows.  One,  computer  power  has  reached  the  point  w  here  a  raidlatitude  model  can 
be  run  globally  for  many  cases  without  using  an  inordinate  amount  of  computer  resources. 
Two,  as  automatic  processing  of  digisonde  data  has  become  available  it  is  becoming 
feasible  to  assemble  global  data  sets  of  fof2  and  hmf2.  While  global  data  sets  of  fof2  have 
been  available  for  years,  the  labor  intensive  processing  required  to  make  true  height 
analysis  made  it  prohibitive  to  assemble  such  data  sets  of  hmf2.  Having  global  datasets  in 
both  quantities  will  serve  as  an  excellent  test  of  how  well  any  model  can  reproduce  an 
actual  midlatitude  ionosphere. 


8.  MODELING  THE  HIGH-LATI TDDE  F-REGION 

In  the  high  latitude  F-region,  we  were  interested  in  developing  a  model  to  study  the 
production  and  evolution  of  polar  cap  patches  and  auroral  blobs.  With  that  in  mind,  we 
began  updating  an  existing  high  latitude  model. 
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S.l  Model  Modifications 


The  first  moditicLiOi  >iis  uivoised  ^  ilu;!  v.c  t!:c\  iuul  .cipiict!  to  uui 

inidhititiid?  tiKxicl.  'Ities  inclufied  updau'd  tiulLisisin  ‘iL-isu.’  Kites,  .c. 

improved  photoionization  source,  a  ntgiiitinie  'oan  e  of  iofn  ■  itnci,  ..m  i  i  ncvv  ssnid  nivsici 
(ltWV187).  Besides  ujidating  various  niicrosCi.?p:c  pararKci-  a.o.u  imth  seveisi!  oi  mi 
inputs,  we  also  conipieteh  rewrote  the  outer  lasers  of  iiie  'OwKii.  ’.)u:'  j-m;'  sc  ssas  to 
develop  a  comprehensive  and  Hcxiblc  time-dcpeedciH  miHsel  di  n  v-rndd  mimv  us  to  assess 
the  roles  of  various  pnKcsses  tn  creating  large-scale  structures  wuliin  the  high  latitude  h 
region.  The  present  nv.x'e!  can  use  any  of  three  higit  iatnude  con*,  c-.  i.iui  nuxieis  i  1  leehs  et 
ah  1982;  Hairston  and  Heeiis,  1990;  and  Heppner  and  Maynard,  i  a-,?)  ii  i,,per;Ucs  in 
what  we  cal!  a  station  mode  svhere  one  selects  a  particular  gro-..nd  f'cam  .n  arid  the  iimc-  ai 
which  you  want  to  calculate  altiuidc  piofiies  of  the  Or  densu;,  diic  proiiiis  arc  caiiuhaeu 
by  following  Bux  tubes  whose  irajeeiories  pass  oscr  the  station  at  tiie  sjvcitied  mnes.  in 
this  way,  you  can  sirmdatc  an  ionosphere  as  it  would  be  (.*bscrvc;i  by  a  gruund-hased 
instrument  such  as  a  digi.sonde.  in  a  first  attempt  to  study  linic  -.  a!  ying  canvection 
patterns,  the  model  can  use  up  to  16  difterent  convection  piutc.ms  during  a  siiigie 
simulation. 


8.2  Theoretical  Studies 

Our  first  study  asing  our  high  laotuJe  inotlcl  consisted  of  cakuhiting  tne  1-2 
parameters  at  Sondrestrom  and  Qanaq  fur  bebruary  i9,  199i)  frttni  K  1  !'  to  just  afte;  17 
L'T.  This  period  was  dictated  by  the  availability  of  IMlh  DMSP,  and  digisonde  data. 
During  this  period,  we  defined  nine  different  convection  patterns  based  on  drift 
measurements  made  by  the  l-H  DMSP  satellite.  Our  first  ealculation  was  for  Sondrestrom 
and  used  ju.st  co-rotaiion  without  any  other  convcctioin  This  was  \o  cdablish  a  btiscline 
with  which  to  assess  the  etfects  of  including  a  conscction  pattern  in  later  eaiculaiions. 
What  we  saw  was  that  K'X'al  production  at  .Sondrc.Mroin  prcxluced  lunpic  ionization  ac 
compared  lo  what  wu'.  observed  but  was  .ety  smooth  compared  m  the  Inghiy  van.ible 
digisonde  data.  We  next  ni.ide  a  calculanon  i.ung  single  convecdcin  n.utcrr.  one  i,f  the 
nine  from  the  pericxl  ol  intec.-n,  and  loiiiuJ  a  i;u.!e  more  temporai  .o.ruwure  as  compaicd  to 
the  C(.>-rotation  results  Ixit  snil  a  much  smoother  resuit  as  compared  to  the  digismuic  data. 
Using  ail  nine  convection  patterns  in  our  simulation  gave  slightly  dif  ferent  results  in  detail 
but  c|ualitatively,  ttiose  results  were  the  .same  .is  tiie  one  patten,  simnl-iiinn 

After  mtiking  a  .similar  >ct  ol' calculations  for  Oanaq  and  obijining  similar  resuus, 
we  made  a  senes  of  calculatu.ms  designed  to  examine  the  sen.oliviiy  of  our  results  to  initial 
conditions,  particle  p.mcipiiation.  and  neutral  wind,  fo  test  scnsiio  ;ty  to  initial  conditions. 
we  repeated  the  nine  convection  pattern  cakuhuion  but  followedi  a  tUix  'uhe's  trajectory  b't 
12  hours  before  it  reached  Soiidrestrom  ratlicr  than  the  'dx  hours  used  in  the  origiita! 
calculation.  This  had  the  eftecf  of  producing  much  mc-re  dramatic  vai  uition  tn  the  temporal 
behavior  in  fof2  at  Sondicstrom.  On  the  other  h.uid,  calculations  indutling  no  particle 
precipitation  showed  very  little  structure  in  eitlior  the  6  or  the  12  hou?  trajectory  cases. 
Turning  off  the  nentrai  wind  likewise  had  an  effect  t.f  .smoothing  out  some  of  ttic  structure. 
.So  what  we  are  seeing  is  that  :n  this  ca.se,  tlie  model  can  produce  .wnne  siructuiv  hut  its 
existence  is  sensitive  to  a  variety  of  parameters  within  the  model.  I'.ut  m  all  cases,  the 
structure  prcxluccd  is  smoother  than  what  is  observed. 

•A  final  set  of  c..dcuiations  was  maiie  that  litcludcd  a  iix-aii/xd  "voricx"  potential 
structure  adder!  to  a  global  pattern  at  a  fixed  location  tn  magricuc  ,ik  admates.  I  he  ct  ieci  n; 
thI^  reetun  of  more  intense  ciecliic  fields  was  to  cause  an  incrc;isc  in  ihc  ('  eifcctisc  iised  ni 


calculating  the  0+  loss  rates.  This  in  turn  leads  to  a  depletion  in  the  F2  peak  density  even 
when  the  vortex  only  existed  for  10  minutes. 


8.3  Future  Directions 

To  make  further  progress,  our  present  model  would  need  two  types  of 
improvements;  (1)  a  more  flexible  and  sophisticated  handling  of  time  varying  convection 
and  particle  precipitation  and  (2)  a  capability  to  do  global  modeling.  When  these 
capabilities  are  available,  then  an  effective  study  can  be  made  of  various  scenarios  on  how 
to  make  F  region  structure. 


9.  Summary 

In  this  study  on  ionospheric  transport,  we  have  examined  a  variety  of  issues  using 
both  kinetic  and  fluid  approaches.  In  our  study  of  electron  flow  between  the  ionosphere 
and  magnetosphere,  we  have  focused  on  two  areas:  a  kinetic  model  for  the  diffu.se  aurora 
and  the  role  of  energetic  photoelectrons  within  the  polar  cap.  The  kinetic  model  of  the 
electron  flux  in  the  diffuse  auroral  region  emphasizes  the  consequences  of  pitch-angle 
scattering,  convection,  and  ionospheric  coupling  on  the  electron  population.  We  found  that 
when  pitch  angle  scattering  is  not  strong,  a  model  with  precipitation  alone  cannot  explain 
the  equatorward  edge  of  the  diffuse  aurora.  In  our  one  dimensional  model,  a  drift- 
boundary  limitation  of  flux  must  be  imposed  in  order  to  produce  a  realistic-looking 
equatorward  edge.  In  our  two  dimensional  model,  a  realistic  edge  can  be  produced  as  the 
results  of  energy-dependent  gradient/curvature  drifts  of  the  electrons. 

In  our  work  on  the  sunlit  polar  cap,  we  have  shown  experimentally  and 
theoretically  that  the  energetic  electron  population  has  an  ionospheric  source 
(photoelectrons)  and  a  magnetospheric  source  (polar  rain).  We  find  that  at  low  energies, 
our  theory  is  in  quantitative  agreement  with  observations  except  for  about  a  factor  of  2  in 
the  low  energy  backscatter.  At  high  energies,  there  is  excellent  quantitative  agreement 
between  theory  and  experiment. 

Our  study  of  ion  flows  between  the  ionosphere  and  magnetosphere  also  focused  on 
two  areas:  a  test  of  the  first  self-consistent  theory  for  the  combined  electron-proton- 
hydrogen  atom  aurora  and  simulations  of  ion  acceleration.  The  test  of  our  coupled  three 
species  auroral  model  consisted  of  calculating  the  upgoing  electrons  that  result  from 
downgoing  protons  incident  at  the  top  of  the  ionosphere.  Measurements  from  the  Low 
Altitude  Plasma  Instrument  on  board  the  Dynamics  Explorer  2  satellite  provided  both  the 
needed  boundary  conditions  and  the  resultant  upgoing  electron  fluxes  with  which  the  model 
results  were  compared.  We  found  that  within  ie  limitations  of  the  data  that  the  agreement 
between  model  and  data  was  good. 

Two  types  of  simulations  were  performed  to  study  ion  acceleration.  First,  two 
dimensional  plasma  simulations  were  used  to  study  the  generation  of  broadband  turbulence 
and  the  resultant  particle  acceleration.  From  these  simulations,  it  was  found  that  a 
broadband  spectrum  of  frequencies,  spanning  the  observed  range  of  auroral  hiss,  could  be 
excited  and  that  ion  acceleration  would  result.  Second,  mesoscale  (.Monte  Carlo) 
simulations  were  used  to  analyze  ion  conic  measurements  from  the  rocket  flight  MARIE. 
What  we  found  was  that  the  observed  turbulence  could  heat  protons  from  1  eV  to  7  eV. 
While  an  underestimation  of  what  was  ob-served,  this  amount  of  heating  is  the  right  order 
of  magnitude. 
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Fluid  models  of  the  E  and  F  region  have  been  developed  to  tonipleineni  our  kinetic 
models.  The  region  model  can  be  used  globally  and  has  been  adjusted  to  produce  a  more 
realistic  FoE.  The  F  region  nuxlel  has  both  a  midlatitude  and  a  high  latitude  version.  The 
midlatitude  version  has  been  tested  against  a  variety  of  measurements  and  has  been  found 
to  reproduce  much  of  the  ionospheric  climatology.  On  the  other  hand,  its  ability  to  specify 
the  F  region  at  a  given  lime  and  place  is  clearly  limited  by  present  uncertainties  asscKiated 
with  the  inputs  to  the  model.  The  high  latitude  F  region  model  is  at  a  much  less  advanced 
stage  of  development  and  there  is  much  to  learn  about  how  it  behaves  under  various 
conditions. 

Over  the  period  of  the  contract,  nine  papers  were  written  and  iwenty-one 
presentations  were  given  at  various  national  and  international  scientific  meetings.  The  nine 
papers  are  listed  below. 

Chang,  T.,  G.B.  Crew.  J.M.  Retterer,  and  J.R.  Jasperse,  Ionic  Conics  and 
Counterstreaming  Electrons  Generated  by  Lower  Hybrid  Waves  in  the 
Earth’s  Magnetosphere,  I  EE  Trans.  Plasma  Sci.,  17,  186,  1989. 

Daniell,  R.E.  Jr.,  D.T.  Decker,  D.N.  Anderson,  J.R.  Jasperse,  J.J.  Sojka. 
and  R.W.  Schunk,  A  Global  Ionospheric  Conductivity  and  Electron  Density 
(ICED)  Model,  in  Proceedings  of  the  6*^  International  Ionospheric  Effects 
Symposium,  p.  351,  1990. 

Decker,  D.T,,  and  J.D.  Winningham,  Energetic  Atmospheric 
Photoelectrons  Due  to  Solar  X-Rays  and  Auger  Production,  in  Physics  of 
Space  Plasmas,  ed.  by  T.  Chang,  G.B.  Crew,  and  J.R.  Jasperse,  SPI 
Conference  Proceedings  and  Reprint  Series,  Vol.  8,  p.  199,  1989. 

Decker,  D.T.,  J.R.  Jasperse,  and  J.D.  Winningham,  Energetic 
Photoelectrons  and  the  Polar  Rain,  in  Physics  of  Space  Plasmas,  ed.  by  T. 

Chang,  G.B.  Crew,  and  J.R.  Jasperse,  SPI  Conference  Proceedings  and 
Reprint  Series,  Vol.  9,  p.  15  1990. 

Reini.sch,  B.W.,  D,  Anderson,  R.R.  Gamache,  X.  Huang,  C.F.  Chen,  and 
D.T.  Decker,  Validating  Ionospheric  Models  with  Measured  Electron 
Density  Profiles,  submitted  to  Proceedings  of  COSPAR  1992,  Wa.shington. 

DC. 

Retterer,  J.M.,  D.T.  Decker,  and  J.R.  Jasperse,  A  Kinetic  Model  for  the 
Diffuse  Aurora,  in  Physics  of  Space  Plasmas,  ed.  by  T.  Chang.  B.  Coppi, 
and  J.R.  Jasperse,  SPI  Conference  Proceedings  and  Reprint  Scries,  Vol.  7, 
p.  167,  1988. 

Retterer,  J.M.,  and  T.  Chang,  Pla.sma  Simulations  of  Intense  VIJ- 
Turbulence  and  Particle  Acceleration  in  the  Suprauroral  Region,  ed.  by  T. 

Chang,  G.B.  Crew,  and  J.R.  Jasperse,  SPI  Conference  Proceedings  and 
Reprint  Series,  Vol.  8,  p.  309,  1989. 

Winningham,  J.D.,  D.T.  Decker,  J.U.  Kozyra,  J.R.  Jasper.se.  and  A.W. 

Nagy,  Energetic  (>  60  eV)  Atmospheric  Photoclectrons,  J.  Gcophvs.  Res., 

94,  15335, 1989. 
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Yasseen,  F.,  J.M,  Retterer,  7’.  Chang,  and  J.D.  Winnmuhani.  Moiue  C'arln 
Modeling  of  Photoelectron  Distributions  with  Aiu>f!uiIous  Heat  Fins. 
Geophys.  Res.  Leu.,  J6,  1023,  1989. 

There  are  five  major  computer  models  that  have  been  developed  under  this  contract 
In  all  five  cases,  they  jire  modified  or  extended  versions  of  earlier  Air  Force  descl-rpcd 
models.  They  are; 

(1)  Photoelectron  transport  model 

(2)  FYoton- hydrogen  atom  transpon  model 

(3)  E  region  model 

(4)  Middle  latitude  F  region  model 

(5)  High  latitude  F  region  model 
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